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ORIGINAL RESEARCH 


Impact of late sowing on morphological and yield traits in 40s bread wheat 


Muhammad Adnan!, Abdullah Khan!, Fida Mohammad!, Fawad Ali*, Quaid Hussain” 


ABSTRACT 
‘Department of Plant Breeding 
and Genetics, The University of The unpredictability and large fluctuation of the climatic conditions in 
rainfed regions influences spring wheat yield and grain quality. These 
variations offer the opportunity for the production of better quality wheat. 
The effect of late sowing on wheat morphology and grain yield was 


Agriculture Peshawar, Pakistan. 
Oil Crops Research Institute, 


reuenea ped oF Beiii studied in different 40s bread wheat at the research farm of PBG, The 
Guac. bianca a a University of Agriculture Peshawar, Pakistan during 2013-14. Forty wheat 


genotypes were tested under normal and late sowing in 5 x 8 alpha lattice 
design with three replicates. Combined analysis of variance exhibited 
significant genotype by environment interactions for days to heading, flag 
leaf area, days to maturity, plant height, spikes m*, grains spike, 1000- 
grain weight, biomass yield, grain yield and harvest index. Days to 
emergence, headings, maturity ranged from 9 to 12, 111 to 121 and 155 to 
164 days under normal while under late planting it ranged from 25 to 29, 


3Department of Plant Sciences, 
Quaid-I-Azam University, 
Islamabad-Pakistan. 
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Mune Enna an 95 tol07 and 137 to 143 days. Mean data under normal planting ranged 
te between 77 to 125m; 25 to 41cm; 99 to 199; 10 to 13 cm 32 to 49; 52 to 
adnan.gxu786@ gmail.com 88g; 8533 to 13667 kg, 1869 to 4681 kg; 21 to 35% whereas under late 
Contact: +8618577966364) planting its range was 63 to 91 cm, 18 to 37 cm’, 57 to 137, 8 to 12 cm, 22 

to 52, 36 to 75g, 2400 to 7933 kg, 540 to 2739 kg and 20 to 42% for plant 
Received: 25 Oct 2021 height, flag leaf area, spikes m~”, spike length, grains spike", 1000-grain 
Revised: 22-Nov-2021 weight, biomass, grain yield and harvest index, respectively. Wheat 
Accepted: 15-Jan-2022 genotypes planted at late condition took maximum days to emergence, 


while less number of days were reacquired for wheat genotypes planted at 
normal sowing date to get mature. Late planting negatively affected all 
yield contributing traits like; spikes m? (29%), grains spike! (18%) 1000- 
grain weight (29%), biomass (55%) and grain yield (50%). On the basis of 
the current exploration, it is obtained that genotype SRN 19111 was 
identified superior for 1000-grain weight, biomass yield and grain yield 
under normal planting, while genotype PR-107 exhibited higher grain 
yield under late planting. Therefore, these genotypes are recommended for 
further extensive testing. 


KEYWORDS: Alpha lattice design, genotypes, sowing time, wheat yield 
and yield components 


This is an open-access article published by the Journal of Soil, Plant and Environment, which permits use, 
distribution and reproduction in any medium, provided the original work is properly cited. 


1. INTRODUCTION its high nutrient values (Adnan et al., 2017; 
Dutta et al., 2019). Cereal crops; wheat, rice 
and maize account for roughly half of the 
world’s human caloric intake (Rosenthal and 


Wheat (Triticum aestivum L.) is one of the 
major cereal crops in the world known for 


www.jspae.com 1 


Journal of soil, plant and Environment 


Ort, 2012). Wheat is regarded as a staple 
food crop in most of the countries including 
Pakistan. World wheat production during 
2009-10 decreased to 0.32% as compared to 
the previous cropping season (2008-09), 
whereas 5.4% decrease in production was 
estimated for the 2010-2011 season. This 
decrease in wheat production might be due 
to several reasons such as_ improper 
agronomic practices, poor management and 
unfavorable weather conditions such as high 
temperature, drought and salinity (Ali et al., 
2019). The assessment and characterization 
of plant germplasm are vital for usage of 
genetic resources of any type (Ali et al., 
2019; Vigna et al., 2011). Wheat genotype 
developed through modern technologies is 
having good yield and achieved to having 
good harvest index. However, there is still a 
gap to develop and select new wheat 
genotype which can resist the harsh 
environmental condition and giving high 
yield (Riazuddin et al., 2010). It has been 
observed that about 80 % of the wheat crop 
in Pakistan is late sown, while only 20 % of 
the wheat crop in sown at normal planting 
time. Planting of the wheat crop at normal 
sowing time could add about two million ton 
of wheat to the national production of wheat. 
The causes of delay in sowing of wheat may 
be due to the presence of previous crops in 
the field and the absence of quality and 
timely agriculture input because of these 
reasons at grain filling stage the wheat crop 
is induced to terminal high temperature due 
to which the yield of the wheat crop is 
reduced much. Production of tillers, kernel 
size, 1000-grain kernel weight, spike length, 
biomass yield, grain yield and harvest index 
are decreased that may vary from genotype 
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to genotype. High temperature may also 
negatively affect other biochemical and 
physiological functions in wheat crop 
(Hamam and Khaled, 2009; Kattenberg et 
al., 1995; Reynolds et al., 2001). Delay in 
sowing of wheat after 10 November may 
cause 42 kg ha'loss in yield (Khan, 2003). 
The probability of invasion of pest disease 
attacks, drought and _ high-temperature 
shocks are more for wheat genotype sown 
late. Therefore, the normal planting of wheat 
crop shows good results for obtaining good 
yield. Molecular and conventional 
approaches can play a vital role in the 
development of wheat genotypes that are 
high-yielding and having adoption to 
various environmental risks. 


Production considerations like vegetative 
growth, grain yield and quality is more 
affected by planting time. The research work 
carried out on planting date may be helpful 
for wheat growers to obtain optimum yield 
from their farming. Late sowing increases 
the risk of yield loss (Ehdaie et al., 2001). 
Morphological and yield traits are directly 
related with the sowing time. Mid- 
November sowing produces the highest 
number of tillers m*, spike m7, 1000 grain 
weight and grain yield (Nasser 2009). Early 
or late sowing increases the risk of yield 
losses (Ehdate et al., 2001). Similarly 
biomass accumulation, grain yield, number 
of spikes m! and 1000 grain weight of 
wheat were increased with early (early 
November) sowing over late (December) 
sowing as reported by Aftab et al. (2004). 
The present study was carried out to find out 
the impacts of late sowing in morphological 
and yield traits in newly developed wheat 
genotypes. 
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2. MATERIALS AND METHODS 


2.1 Experimental location and weather 
detail 


The field experiment was carried out at the 
Research Farm of PBG, The University of 
Agriculture Peshawar, Khyber 
Pakhtunkhwa, Pakistan during 2013-14. The 
farm is located at 34.01_ N latitude, 71.35_ 
E longitude, at an altitude of 350 m above 
sea level in Peshawar valley. The site is 
situated about 1600 km North of the Indian 
Ocean and has a continental type of climate. 
Warsak canal from Kabal River was used 
for the irrigation system. The soil of this site 
is clay loam, containing low in organic 
matter (0.87%), potassium (121 mg kg’!); 
phosphorus (6.57 mg kg’), alkaline (pH 
8.2); and with calcareous nature (Khan et al. 
2009, So et al. 2012). The minimum, 
maximum temperature and average rainfall 
data was collected from the metrological 
station and shown in Fig | a & b. 
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2.2 Experimental design and agronomic 
practices 

Experimental material consisted of 40 
genotypes (39 Advance lines and 1 local 
check). These genotypes were planted under 
normal and late planting dates to study the 
effects of planting time on _ different 
genotypes. Experimental material was 
received from the National Agriculture 
Research Center (NARC) Islamabad as 
National Uniform Yield Trial (NUYT) 
2013-14. The experiment was laid out in 5 x 
8 alpha lattice design with 3 replicates. 
Every genotype was planted in a 6-row plot, 
having row to row and plant to plant 
distance of 0.25 m and 5 m, respectively. 
Normal planting was done on November 20, 
2013, while late planting was done on 
December 21, 2013. All other cultural 
practices were applied uniformly to the 
experiments sown under normal and _ late 
planting dates. Three to four irrigations were 
applied during the wheat growth period. 
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Figure 1. Mean monthly precipitation and air temperature during Nov 2013 to August 2014. 
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Table. 1. List of genotypes used in the experiment 


Genotype Code Breeding center 
1 109384 RARI- Bahawalpur 
2 99172 RARI — Bahawalpur 
3 99346 RARI — Bahawalpur 
4 99114 RARI — Bahawalpur 
5 DN-93 ARI- DI Khan 
6 CT 09137 NIFA-Peshawar 
i) SRN 09111 NIFA — Peshawar 
8 V-09082 WRI — Faisalabad 
9 V-09087 WRI — Faisalabad 
10 V-10104 WRI — Faisalabad 
11 V-10110 WRI — Faisalabad 
ibe. V-11160 WRI — Faisalabad 
13 SKD-11 WRI-Sakrand 
14 NN-Gandum-I NIBGE Fsd 
15 NN-Gandum-II NIBGE Fsd 
16 TW96010 AZRI Bhakkar 
17 TW96018 AZRI Bhakkar 
18 SD-998 NIA Tandojam 
19 NIA-MN-08 NIA Tandojam 
20 CIM-04-10 NIA Tandojam 
21 ESW-9525 NIA Tandojam 
22 PR-103 CCRI-Pirsabak 
23 PR-106 CCRI-Pirsabak 
24 PR-107 CCRI-Pirsabak 
25 RCA-1 RCA Seeds 
26 V-11005 WRS-Tandojam 
27 NR-413 NARC-Islamabad 
28 NR-421 NARC-Islamabad 
29 NR-409 NARC- Wheat 
30 NR-419 NARC-Islamabad 
31 UAF-9452 Univ. of Agri. Fiasalabad 
32 Guard-C Hybrid — Guard 
33 SAWSN-02-102 AZRC-DI Khan 
34 Janbazz UOA Peshawar 
33 TD-1 WRI-Sakrand 
36 Pirsabak-13 CCRI-Pirsabak 
37 Sehar-06 WRI-Faisalabad 
38 V07096 WRI-Faisalabad 
39 Aas-11 RARI- Bahawalpur 
40 NARC-11 Wheat-NARC 
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2.3 Data collection and measurements 


Days to emergence of wheat were 
observed in all plots by counting the number 
of days from sowing until the date when 
more than 70% of plants emerged. Days to 
heading were recorded by counting days in 
number from sowing until the date when 
more than 75% of tiller formed anthers. Data 
on plant height was recorded by taking the 
height of 5 tillers at random in each plot and 
the height from the base to the tip of each 
tiller was measured using a meter rod 
including awns and was then averaged. Flag 
leaf area was recorded at the panicle 
initiation stage. Days to physiological 
maturity of wheat were observed in all plots 
by counting the number of days from 
sowing until the date when more than 75% 
of tillers had become mature in all plots. 
Spikes m7? was recorded by counting spikes 
in three middle rows and were then 
converted to m®” using the same formula as 
spikes m*. Spike length was recorded at 
physiological maturity by using simple 
geometric ruler. Data on leaf chlorophyll 
content was recorded as a SPAD value using 
the Spad meter by putting the flag leaf on 
the scanner of the meter and hold for a while 
(Islam et al. 2014). Grains spike’! was noted 
by counting number of grains per spike in 
five randomly selected spikes in each plot 
and was averaged. Data on thousand grain 
weight was obtained by counting 1000 
grains at random from the grain lot of each 
plot and were measured using an electronic 
balance. Biological yield was noted by 
harvesting three rows in the middle of each 
plot and were dried in an open field for one 
week and were then weighed and converted 
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to kg ha"! using equation 1. The grain yield 
was determined by threshing the harvested 
sample and the grain obtained was weighed 
and converted to kg ha’!. Harvest index 
expressed in percentage is the ratio of seed 
yield to biological yield. It was determined 
by dividing seed yield by biological yield 
multiplied by 100. 

BY(kg ha‘!)= BY(kg) in four central rowsx 
10000/rows-row distance xX number of 
number of rows X row length------------- (1) 


2.3 Statistical Analysis 

Analyses were carried and the significant 
means for various traits were separated with 
the application of LSD test. Sigma plot X7 
was used for figure and data analysis was 
done through IBM-SPSS20 and Microsoft 
excel. 


3. RESULTS 
3.1 Phenology 


Data on days to emergence, days to 
heading and days to maturity was 
considerably altered by wheat genotypes, 
normal and late planting. Days to emergence 
for wheat genotypes was recorded from 9 to 
12 days under normal planting sowing and 
25 to 29 days under late planting sowing 
condition (Table 2). Whereas among 
genotypes 109384, DN-93, V-10104, V- 
10110, V-11160 and NN-Gandum-I took 
minimum (9) days to emergence while 
genotype SD-998 took maximum (12) days 
to emergence. In contrast, genotypes DN-93 
and UAF-9452 took minimum (25) days to 
emergence while genotype ESW-9525 took 
maximum (29) days to emergence under late 
planting condition. Average over normal and 
late planting days to emergence ranged from 
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-® Treatment vs Early sowing 
-O- Treatment vs Late sowing 


Chlorophyll (SPAD) 


INA 


My Nia 


Genotypes 


Figure 2. Chlorophyll content as influenced in different wheat genotypes and by sowing dates. Note: the serial 


numbers from | to 40 indicates the genotypes, for detail see table 1. 


18 days for genotype NR-419 to 26 days for 
genotype PR-107. Mean values for 40 wheat 
genotypes for days to emergence were 10 
and 26 days under normal and late planting 
sowing conditions, respectively (Fig. 3A). 


Data on days to heading showed 
minimum (111) days to heading were 
observed for genotype PR-419 and 
maximum (121) days to heading were 
observed for genotypes CT 09137, SD-998 
and NIA-MNO8 under normal planting 
while at late planting minimum (95) days for 
genotypes TW96018 and PR-419 and 
maximum (107) days to heading were 
observed for genotypes DN-93 and ESW- 
9525 (Table 2). Genotype PR-419 took 
minimum days to heading under both 
normal and late planting conditions. The 
significant GxE interaction implies that the 
genotypes were observed to have different 
relative days to heading across planting 
dates. Mean genotypic performance across 
normal and late plantings for days to 
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heading were 118 days in normal planting 
while 103 days in late planting(Fig. 3B). 
Thus in late sowing conditions the head in 
emergence 15 days were reduced as 
compared to normal sowing. 


Data on day’s physiological maturity 
under normal planting ranged from 155 to 
164 days, while under late planting 
condition it ranged from 137 to 143 days 
(Table 2). Genotype UAF-9452 took 
minimum (155) days to physiological 
maturity while maximum (164) days to 
maturity were recorded for genotype NR- 
413 under normal planting. Under late 
planting condition minimum (137) days to 
maturity were recorded for genotype UAF- 
9452 while maximum (143) days to maturity 
were recorded for genotype NR-413. 
Interestingly the genotype UAF-9452 took 
minimum days while genotype NR-413 took 
maximum days to maturity under both 
normal and late planting. Mean days to 
maturity across two production 
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environments ranged from 147 for genotype 
99346 to 153 for genotype V-09082. 
Average days to maturity were 159 and 140 
days under normal and late planting 
conditions, respectively (Fig. 3C). From the 
net difference in days to maturity it was 
notice that the wheat genotypes planted at 
late planting condition get mature 19 days 
earlier than wheat genotypes sown at their 
normal sowing condition. 


3.2 Growth traits 


Data on wheat plant height were 
recorded significantly among genotypes and 
planting dates (Table 2). Generally, shorter 
plants were observed in late sowing as 
compared to the normal sowing date. Under 
normal planting conditions minimum (77 
cm) plant height was observed for genotype 
TD-1 while genotype SAWN-02-102 
showed maximum (125 cm) plant height. 
While genotype TD-1 showed minimum (63 
cm) plant height and genotype SAWN-02- 
102 maximum (91 cm) plant height under 
late planting. Plant height in overall 
interaction value ranged from 63 cm for 
genotype TD-1 to 125 cm for genotype 
SAWN-02-102. Mean values across normal 
and late plantings were observed for 
genotype TD-1 to be 70 cm and 108 cm for 
genotype SAWN-04-102. Under normal 
planting the mean value was 105 cm while 
under late planting mean plant height of 80 
cm was observed (Fig. 3D). 


Spike length under normal planting varies 
from 10 to 13 cm and 8 to 12 cm under late 
planting condition (Table 2). Under normal 
planting genotype Guard-C showed 
mininum (10 cm) spike length while 
genotype 109384 showed maximum (13 cm) 
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spike length. In contrast, genotype CT 
09137 showed minimum (8 cm) spike length 
under late planting condition while genotype 
Aas-11 showed maximum (12 cm) spike 
length. Therefore, in normal sowing the 
genotype produced plants with maximum 
spike length. Mean value for spike length in 
overall interaction was ranged from 8 cm for 
genotype CT 09137 to 13 cm for genotype 
109384. Mean over the two_ planting 
environments genotype Aas-11 produced 
plants with maximum spike length 13 cm 
while genotype CT 09137 produced plants 
with minimum spike length 9 cm. Mean for 
40 wheat genotypes under normal and late 
planting were 12 and 10 cm respectively 
(Fig. 3E). 


Days to Emergence Daysto Heading 


| | 


Days to Maturity Plant Height 


w C “D 


Flag Leaf Area (cm?) Tiller (m7?) 


WE |e 


Figure 3. Impact of normal and late planting 
on yield and yield components of 40 breads 
wheat genotype. 


Maximum (199) number of spikes m7 
were recorded for normal planting while 
minimum (57) spikes m* for late planting 
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(Fig. 3F). Genotype Guard-C showed 
maximum (199) spikes m? while genotype 
SAWSN-02-102 showed minimum (99) 
spikes m* under normal planting condition. 
Under late planting condition genotype 
CIM-04-10 showed minimum (57) number 
of spikes m? whereas genotype PR-107 
showed maximum (145) number of spikes 
m”. The interaction mean value for number 
of spikes m7 for normal and late planting 
ranged from 166 for genotypes Guard-C to 
81 for genotype CIM-04-10 (Table 2). Mean 
values for 40 genotypes for spikes m” over 
normal and late planting were 140 and 100 
respectively. 


3.3. Physiological traits 


The data recorded for flag leaf area under 
late planting ranged from 18 cm? to 37 cm’, 
while under normal planting condition it 
ranged from 41 to 25 cm* (Table 2). 
Genotype DN-93 showed minimum flag leaf 
area of 18 cm’, while genotype V-11005 
exhibited maximum flag leaf area of 37 cm? 
under late planting while in normal planting 
condition minimum (25 cm’) flag leaf area 
were observed for genotype NR-419, while 
maximum (41 cm’) flag leaf area were 
recorded for genotypes PR-107 and RCA-1. 
The significant G x E interaction can be 
confirmed from the mean values for flag leaf 
area as the genotype showing maximum flag 
leaf area in normal sowing did not show 
maximum flag leaf area in late sowing 
which meant that the ranking of genotypes 
changed on the basis for their performance 
for flag leaf area under normal and _ late 
plantings. Based on the mean performance 
of genotypes under normal and late planting 
minimum flag leaf area was observed for 
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genotype V-11160 (25 cm’) while maximum 
for genotypes PR-107 and RCA-1 (38 cm’). 
Averaged over 40 genotypes, flag leaf area 
was 27 and 34 cm’ under normal and late 
planting conditions, respectively. 


Data on chlorophyll content significant 
difference between genotypes and sowing 
dates. Between sowings dates early sowing 
showed higher chlorophyll content across 
the genotypes (Fig. 2). Among the 
genotypes no 7 in the Table 2 (SRN 09111) 
showed higher (51.1) chlorophyll content at 
early sowing (Figure 2), whereas during late 
sowing genotype number 24 resulted higher 
(50.6) chlorophyll content. The possible 
reason for this variation might be attributed 
to the genetic characteristics of the genotype 
as reported by (Kochak-Zadeh et al. 2013) 


3.4 Yield components and yield traits 


Data on number of grains spike’! were 
found significant for genotypes and planting 
date. Minimum number of grains spike’! 
(32) were observed in late sowing as 
average across 40 wheat genotype while 
normal sowing showed maximum (39) 
number of grains spike! (Fig. 3F). Under 
late planting condition minimum (22) 
number of grains spike’ was observed for 
genotype CT 09137, while genotype Janbaz 
showed maximum (52) number of grains 
spike'. In contrast genotype CIM-04-10 
produce maximum (59) number of grains 
spike while genotype Guard-C produce 
minimum (27) numbers of grain spike under 
normal planting (Table 2). Planting dates 
showed significant variation in number of 
grains spike’! in wheat genotypes. Number 
of grains spike’ in overall interaction value 
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ranged from (22) for genotype CT 09137 to 
(59) for genotype CIM-04-10. Due to late 
planting 18% reduction occurred in number 
of grains spike”. Sial et al. (2005) reported 
reduction in grain numbers due to late 
sowing and temperature Stress. The 
significant G x E interaction can be 
confirmed from the mean values for number 
of grains spike! as the genotype showing 
maximum number of grains spike! in 
normal sowing did not produce plants with 
maximum number of grains spike” in late 
sowing which means that the ranking of 
genotype in changed on the basis for their 
performance for number of grains spike"! 
under normal and late planting. 


Minimum (8533) biomass was recorded 
for genotype TD-1 and maximum (13667) 
biomass for genotype SRN 09111 at normal 
planting condition while at late planting 
minimum (2400 kg) and maximum (7933kg) 
biomass for genotype SD-998 and NR-421 
respectively. Maximum and minimum 
biomass was observed for normal and late 
plantings which showed that the interaction 
effect of genotype and normal and _ late 
planting was effective in causing variation in 
genotypes for biomass. The genotypic effect 
was also observed significant which show 
that the genotypes are significantly different 
from each other for biomass. Averaged over 
normal and late planting minimum (5767 kg 
ha!) biomass was observed for genotype 
TD-1 while maximum (9700 kg ha") for 
genotype NR-421. At late planting condition 
55% reduction in bio-mass_ yield was 
observed. The 40 wheat genotypes produce 
an average 10885 kg ha! at normal 
plantings while the same genotypes sown 
late produce an average yield of 5990 kg ha” 
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' (Fig. 4K). It is obvious from the result that 
the wheat genotypes sown at their normal 
sowing period produce biomass yield of 
4895 kg ha! than the wheat genotypes sown 
at late condition. 


Data on 1000 grain weight demonstrates 
slighter (52 g) weight for genotypes SD-998 
and NIA-MN-08, whereas heavier (88 g) 
1000 grains weight was recorded for 
genotype SRN 09111 at normal planting 
while at late planting (36 g) and (75 g) for 
genotypes NIA-MN-08 and V_ 07096 
respectively (Table 2). Interestingly the 
genotype NIA-MN-08 _ has _ resulted 
minimum 1000-grain weight under normal 
and late plantings. In overall interaction 
value for 1000 grain-weight ranged from (36 
g) for genotype NIA-MN-08 to (88 g) for 
genotype SRNO9111. The mean calculated 
for normal and late planting for 1000-grain 
weight genotype NIA-MN-08 had minimum 
(44 g) 1000-grain weight of and genotype V 
07096 had maximum (75 g) of 1000 grain- 
weight. The average value calculated for 40 
wheat genotype was 70 and 50g for normal 
and late planting respectively (Fig. 4G). 
Thus the kernel weight of wheat genotypes 
was reduced grain because of late sowing 
condition and hence normal planting was 
superior to late planting. The averaged 
1000-grain weight of wheat genotypes was 
20 g more than the wheat genotypes sown at 
late planting condition. 29 % weight loss 
was noticed at late planting condition. 


Data on grain yield ranged averaged from 
540 kg ha! for genotype NIA-MN-08 to 
5336 for genotype SRN 09111 (Table 2). 
Data showed that under normal planting 
genotype NIA-MN-08 resulted Minimum 
(1869 kg ha’) grain yield and genotype 
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SRN 09111 resulted in higher grain yield, 
while at late planting genotype NIA-MN-08 
demonstrate lesser grain yield (540 kg ha”) 
and genotype PR-107 showed greater (2739 
kg ha’') grain yield. Across normal and late 
plantings smallest amount of (1205 kg ha‘') 
grain yield kg ha! was confirmed for 
genotype NIA-MN-08 and highest amount 
of (3329 kg ha’) grain yield for genotype V- 
11005. Normal planting produce more grain 
yield kg ha’! than late planting (Fig. 4K). 


1000 grains weight (g) Grain Yield (kgha*) 


"G “HH 


Spike Length (cm) Biomass (kg ha) 


7 my 


Grain Yield (kgha*) Harvest Index (%) 


| | : 


Figure 4. Impact of normal and late planting 
on yield and yield components of 40 breads 
wheat genotype. 

Harvest index under normal planting 
ranged from 21 to 35 percent, while under 
late planting the harvest index ranged from 
20 to 42 percent (Table 2). Genotype RCA-1 
gave maximum harvest index of 35 percent 
whereas genotype SD-998 shows minimum 
harvest index of 21 percent under normal 
planting. Under late planting condition 
genotypes SD-998 and Guard-C_ gave 
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maximum harvest index of 42 percent while 
genotype gave minimum harvest index of 20 
percent (Fig. 4L). 


4. DISCUSSION 

Planting time is one of the major causes of 
low yield in wheat. In normal course, the 
sowing of wheat must be completed in the 
month of November in Pakistan. But in 
some cases, wheat grown after Cotton, Rice, 
Sugarcane and fodder crops in rotation may 
be delayed up-to December. This reveals 
negative impact on the grain yield in wheat. 
Research activities regarding the role of 
planting time and identification of wheat 
genotypes having better performance in such 
situation is highly desirable. Whereas among 
genotypes 109384, DN-93, V-10104, V- 
10110, V-11160 and NN-Gandum-I took 
minimum (9) days to emergence while 
genotype SD-998 took maximum (12) days 
to emergence. In contrast, genotypes DN-93 
and UAF-9452 took minimum (25) days to 
emergence while genotype ESW-9525 took 
maximum (29) days to emergence under late 
planting condition. Average over normal and 
late planting days to emergence ranged from 
18 days for genotype NR-419 to 26 days for 
genotype PR-107. Mean values for 40 wheat 
genotype for days to emergence were 10 and 
26 days under normal and late planting 
sowing conditions, respectively. Whereas 
genotype PR-419 showed less (111) days to 
heading(111) and more (121) days to 
heading were recorded for genotypes CT 
09137, SD-998 and NIA-MNO8_ under 
normal planting, while at late planting 
minimum (95) days for genotypes TW96018 
and PR-419 and maximum (107) days to 
heading were observed for genotypes DN-93 
and ESW-9525 (Table 2). 
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Table 2: Mean performance of wheat genotypes for spikes m7, spike length, grain spike”, 1000 
grain weight, bio-mass yield kg ha", grain yield kg ha"! and harvest index for wheat genotypes 


under Normal and late plantings 2013-14. 


Normal Late 
Traits Ranges | Means | Best LSD Ranges | Means | Best LSD 
Genotype | (0.05) Genotype | (0.05) 
Days to 9-12 10 109384 0.8 25-29 | 26 UAF- 1 
emergence 9452 
(no) 
Days to 95-121 | 118 NR-409 | 0.8 95-107 | 103 NR-409 | 1.4 
heading 
(no) 
Plant 76.93- | 104.41 | V-09082 | 1.2 63.26- | 80.38 Janbaz 3.4 
Height (cm) | 125.2 90.66 
Flag leaf 25.26- | 34.40 | PR-107 1.6 17.69- | 27.28 V-11005 | 2.1 
Area (cm’) | 40.86 36.73 
Days 155- 159 UAF- 4.6 137- 140 UAF- 1 
maturity 164 9452 143 9452 
(no) 
Spike m7 99-199 | 140 Guard-C | 7 56-145 | 100 CIM-04- | 7 
(no) 10 
Spike 9.7- 11.63 | 109384 0.7 8.31- 9.92 Aas-11 0.6 
length (cm) | 13.2 11.56 
Grains 27-59 | 39 CM-04- | 2 22-53 | 32 Janbaz 2 
spike’! 10 
1000- grain 51.74- | 70 SRN 1.9 36.27- | 50.24 V07096 | 2 
weight (g) 87.59 09111 75.46 
Riguiiase 8533- | 11089 | SRN 403 2400- | 4895 NR-421 | 390 
yield (ke 13667 09111 7933.3 
ha’) 
eae 1869- | 3127 SRN 51 540- 1548.21 | PR-107 68 
ae 5336 09111 27389 
Harvest 21-35 | 28 RCA-1 1.4 20.4- 32.4 V-9082 pe | 
index 41.6 


While at late planting minimum (95) days 
for genotypes TW96018 and PR-419 and 
maximum (107) days to heading were 
observed for genotypes DN-93 and ESW- 
9525 (Table 2). Genotype PR-419 took 


minimum days to heading under both 
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normal and late planting conditions. Data on 
days to physiological maturity under normal 
planting ranged from 155 to 164 days, while 
under late planting condition it ranged from 
137 to 143 days. Genotype UAF-9452 took 


minimum (155) days to physiological 
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maturity while maximum (164) days to 
maturity were recorded for genotype NR- 
413 under normal planting. Under late 
planting condition minimum (137) days to 
maturity were recorded for genotype UAF- 
9452 while maximum (143) days to maturity 
were recorded for genotype NR-413. 
Average days to maturity were 159 and 140 
days under normal and _ late planting 
conditions, respectively. It was observed 
that wheat genotypes sown at late planting 
condition emerged 16 days later than wheat 
genotypes sown at normal planting time to 
normal planting condition. The reason 
behind maximum days to emergence may be 
due to low temperature at late planting due 
to which wheat genotypes took maximum 
days to emerge as compared The same 
report was given by Benjamin (1990) and 
Gul et al (2012) as they observed that low 
temperature during emergence and seedling 
growth has detrimental effect on the crop 
establishment and _ productivity. Early 
emergence of head in wheat is one of the 
prime objectives in breeding programs of 
wheat because maximum time is available 
for grain filling. Late heading results in 
decrease grain size resulting in low grain 
weight (Irfaq et al, 2005). Among wheat 
genotype highly significant differences 
between genotype and _— environment 
interaction for days to heading are similar 
with the results of Razzaq et al. (1986), 
Subhan et al. (1991), Inamullah et al. 
(2007), Ilyas et al. (2013). Muhammad et al., 
(2007) conducted a research and observed 
significant result for days to heading across 
two environmental conditions. From the net 
difference in days to maturity it was notice 
that the wheat genotypes planted at late 
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planting condition get mature 19 days earlier 
than wheat genotypes sown at their normal 
sowing condition. The wheat genotype took 
minimum days to get mature at late sowing 
because of continuity of race and to avoid 
from upcoming unfavorable condition. 
Nahar et al. (2010) reported up to 15% 
reduction in maturity period of wheat 
genotypes due to the effect of heat stress. 
Ilyas et al (2013) also found significant 
differences among wheat genotypes for days 
to maturity. Muhammad et al. (2007) also 
reported highly significant genotype by 
planting date interaction for days to 
maturity. The wheat genotype took 
minimum days to get mature at late sowing 
because of continuity of race and to avoid 
from upcoming unfavorable condition. 


4.2 Growth traits 


Plant height, Spike length and number of 
spikes m* were recorded significantly 
among genotypes and planting dates. 
Generally, shorter plants were observed in 
late sowing over normal sowing date. 
Genotype TD-1 results in shorter (77 cm) 
and genotype SAWN-02-102 showed taller 
(125 cm) plant height plant height under 
normal planting condition. Whereas under 
late planting genotype TD-1 showed 
minimum (63 cm) plant height and genotype 
SAWN-02-102 maximum (91 cm) plant 
height. Spike length under normal planting 
varies from 10 to 13 cm and 8 to 12 cm 
under late planting condition. Under normal 
planting genotype Guard-C showed 
minimum (10 cm) spike length while 
genotype 109384 showed maximum (13 cm) 
spike length. In contrast, genotype CT 
09137 showed minimum (8 cm) spike length 
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under late planting condition while genotype 
Aas-11 showed maximum (12 cm) spike 
length. Number of spikes m* were recorded 
from 57 to 199 spikes m*?. Among 
genotypes, genotype Guard-C showed 
maximum (199) spikes m? while genotype 
SAWSN-02-102 showed minimum (99) 
spikes m? under normal planting condition. 
The possible reason for increasing in plant 
height among cultivars might be due to the 
differences in their genetic makeup. These 
results are in agreement with those obtained 
by Wahid et al. (2017). Genotypic 
differences were also observed by Ahmad et 
al. (1997). Similarly, Laghari et al. (2012) 
found similar result that 32.54% reduction in 
plant height took place due to late sowing 
when compares to normal sowing. 
Reduction of 25 cm occurred as a result of 
late planting. Rashid et al. (2004) and Knapp 
and Knapp (1978) have also observed that 
wheat genotypes sown at late planting 
results in reduction in plant height. Irfaq et 
al., (2005) also reported reduction in plant 
height of wheat genotypes due to late 
sowing and high temperature stress. 


Difference in spike lengths of the 
genotypes was also reported by Kakar et.al. 
(2003). Gul et al (2012) also found that 
sowing of wheat genotype at normal time 
produce maximum number spikes m” than 
late sowing. These significant differences 
are pointing towards the presence of 
sufficient genetic variations among the 
genotype. The results are in accordance with 
the findings of Munir et al (1999) and 
Rajora (1999). 


Flag leaf area under late planting ranged 
from 18 cm? to 37 cm’, while under normal 
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planting condition it ranged from 41 to 25 
cm’. Averaged over 40 genotypes, flag leaf 
area was 27 and 34 cm? under normal and 
late planting conditions, respectively 
.Genotype DN-93 showed minimum flag 
leaf area of 18 cm’, while genotype V-11005 
exhibited maximum flag leaf area of 37 cm? 
under late planting while in normal planting 
condition minimum (25 cm”) flag leaf area 
were observed for genotype NR-419, while 
maximum (41 cm’) flag leaf area were 
recorded for genotypes PR-107 and RCA-1. 
Similarly for chlorophyll content, among the 
genotypes SRN 09111 showed higher (51.1) 
chlorophyll content at early sowing (Figure 
2), whereas during late sowing genotype 
PR-107 resulted higher (50.6) chlorophyll 
content. The possible reason for this 
variation might be attributed to the genetic 
characteristics of the genotype as reported 
by Kochak-Zadeh et al. (2013). Rane et al. 
(2007) reported that cooler climates favored 
the vegetative as well as the reproductive 
phases of wheat growth. 


4.4 Yield components and yield traits 


In wheat, grains per spike, biomass, 1000 
grain weight and grain yield is one of the 
important primary yield contributing 
character. Number of grains spike”, 
biomass, 1000 grain weight and grain yield 
were found significant for genotypes and 
planting date. Under late planting condition 
minimum (22) number of grains spike! was 
observed for genotype CT 09137, while 
genotype Janbaz showed maximum (52) 
number of grains spike’. Similarly less 
biomass was recorded for genotype TD-1 
and higher biomass for genotype SRN 
09111 at normal planting condition while at 
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late planting minimum and maximum 
biomass were recorded for genotype SD-998 
and NR-421 respectively. The 40 wheat 
genotypes produce an average 10885 kg ha’! 
at normal plantings while the same 
genotypes sown late produce an average 
yield of 5990 kg ha‘. Data on 1000 grain 
weight demonstrates slighter weight for 
genotypes SD-998 and NIA-MN-08, 
whereas heavier 1000 grains weight was 
recorded for genotype SRN 09111 at normal 
planting while at late planting lower and 
higher 1000 grains weight for genotypes 
NIA-MN-08 and V_ 07096 respectively. 
Grain yield showed that under normal 
planting genotype NIA-MN-08 resulted 
minimum grain yield and genotype SRN 
09111 resulted in higher grain yield, while at 
late planting genotype NIA-MN-08 
demonstrate lesser grain yield and genotype 
PR-107 showed greater grain yield. Across 
normal and late plantings smallest amount of 
grain yield kg ha! was confirmed for 
genotype NIA-MN-08 and highest amount 
of grain yield for genotype V-11005. The 
possible reason in the differences of bio- 
mass and yield attributes of the genotypes 
might be attributed to the genetic 
characteristic of each genotype as reported 
by Kakar et al. (2003). Akram et al. (2008) 
also found significant variation in wheat 
genotypes for 1000-grain weight. Abdullah 
et al. (2007) and Ansari et al. (1989) also 
found that late planting of wheat genotypes 
results in decrease of 1000- grain weight. 
Laghari et al (2012) and Irfaq et al (2005) 
reported biomass yield reduction in wheat 
genotypes because of late sowing. The 
finding of our research that the ranking of 
wheat genotypes according to _ yield 


www.jspae.com 


Adnan et al. 


production are changed by sowing the same 
wheat genotypes at two _ different 
environment are similar with the research 
findings of Cotes et al. (2006), Amin et al. 
(2005), Khalil et al. (2005), Garcya et al. 
(2003) and Kakar et al. (2003). Rane et al. 
(2007) reported that colder climates favored 
the vegetative as well as the reproductive 
phases of wheat growth 


5. CONCLUSION 


From this study it is concluded that late 
sowing of wheat negatively affect both 
morphological traits and grain yield and may 
results in 50 % of wheat yield loss. Hence 
timely sowing of wheat genotypes is highly 
suggested. The conducted research 
identified genotype SRN 19111 superior for 
1000-grain weight, bio-mass yield and grain 
yield under normal planting, while genotype 
PR-107 exhibited higher grain yield under 
late planting. Therefore, these genotypes are 
recommended for further extensive testing 
and utilization in different wheat breeding 
programs. 
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ABSTRACT 


In Khyber-Pakhtunkhwa, Pakistan, wheat yield is subjected to the 
availability of water and the proper rate of Nano-black carbon in the soil. 
Delay in the rain and unsuitable soil health cause severe yield reduction. 
Therefore this experiment was conducted to compare irrigation levels in 
relation to a different rate of Nano-black carbon to find out the high 
yielding fact that could enhance wheat productivity and food security. 
Three different irrigation levels (250-mm, 275-mm and 300-mm), were 
compared in early growth stages with five different rates of Nano-black 
carbon (5Mg ha'!, 10Mg ha", 15Mg ha", 20Mg ha"! and 25Mg ha’). All 
other agronomic practices were kept similar for each replication. Data was 
recorded on different growth parameters such as days to emergence, 
emergence m”, plant height, spike length, number of spikes m7”, thousand- 
grain weight and grain yield. The study confirmed that almost all irrigation 
levels were prominent but a significant reduction in different parameters 
was observed with variation in Nano-black carbon that could ultimately 
affect soil health and productivity. It was concluded that the proper rate of 
Nano-black carbon can significantly enhance the development of the roots 
system which may ultimately increase the shoot growth and final yield. 
Wheat irrigation levels (250 mm) can properly save water and increase 
wheat productivity in combination with Nano-black carbon amendments. 


KEYWORDS: Wheat (Triticum aestivum L.), Nano-black carbon, 
growth, yield, drought, irrigation regimes 


This is an open-access article published by the Journal of Soil, Plant and Environment, which permits use, 
distribution and reproduction in any medium, provided the original work is properly cited. 


1. INTRODUCTION 


Wheat (Triticum aestivum L.) belongs to 
the family Gramineae and is consumed as a 
major grain crop of the world (Imran, 
2021a). It is also known as the king of cereal 
and supplies about 60% of the calories and 
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protein of the average diet (Khalil and Jan, 
2012). It is ranked 1“ among other cereal 
crops in Pakistan and occupies about 66% of 
the annual food cropped area (Anonymous, 
2012). In Pakistan, wheat is grown on an 
area of 9204 thousand hectares with an 
annual production of 9260 thousand hectares 
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with an average yield of 2752kg ha” (Tunio 
et al. 2016). Despite higher yield potential, 
the average yield of wheat in Pakistan is low 
as compared with advanced countries 
(Imran, 2021a; Imran et al., 2021a; Imran et 
al., 2020a). Several factors are responsible 
for the low productivity of wheat in Pakistan 
like edaphic factors, cultural practices, 
genetic factors and — environmental 
conditions. In the current scenario, improved 
cultural practices and proper nutrient 
management are the major issues in wheat 
production. In Pakistan wheat are raising 
under both, irrigated and rain-fed condition, 
whereas in Khyber Pakhtunkhawa, almost 
70% of wheat is grown under rainfed 
conditiondmran, 2015a; Imran, 2021b; 
Imran et al., 2021c). Efforts are in progress 
to breed for high yield, better protein 
quality, biotic and abiotic stresses resulting 
in the development of many cultivars (AlI- 
Maskri et al., 2012). Irrigation levels having 
the potential for early seedling vigor use 
available soil water efficiently, resulting in 
better dry matter accumulation and higher 
grain yield (Awan et al.,2005). The best 
irrigation level allows the cropto produce a 
higher yield under unfavorable conditions 
(drought, salinity and extreme temperature) 
(Imran et al., 2020b; Imran et al., 2018; 
Imran et al., 2021b). 


Grain yield of wheat can be reduced with 
the reduction in irrigation level, whereas the 
addition of biochar might improve water use 
efficiency and plant essential nutrients 
uptake (Jan et al., 2012). Previously it is 
reported that biochar with different 
irrigation regimes may show high efficiency 
and booted yield (Ali et al., 2020; Ali et al., 
2021; Imran et al., 2020c; Imran et al., 
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2021d; Imran et al., 2021e; Imran, 2015b). 
Donaldson et al. (2013) reported that 
irrigation regimes resulted in increased 
wheat straw production and generally higher 
grain yield compared with various carbon 
sources. Furthermore, Patel et al. (2010) and 
Aslam et al. (2013) tested newly introduced 
irrigation levels, concerning yield and 
observed that the highest yield (5484kg ha‘!) 
was obtained with the integration of organic 
sources and the highest irrigation level 
(Tunio et al. 2016). Likewise, yield 
reduction of 27% and 52% was noted by Ali 
et al. (2014) when wheat crop was sown 
with and without organic amendments. In 
addition, Tahir et al. (2011) concluded that 
regardless of irrigation levels, better yields 
were obtained when wheat was treated with 
organic sources application before irrigation. 
Shah and Akmal (2002) reported that grain 
yield decreased by 28 percent when 
irrigation was delayed by 45 days in the 
season. Furthermore, early irrigation favored 
in high tillering and ultimately in respect of 
grain and biological yield. 


The current study aimed to investigate the 
effects of different irrigation levels with five 
different Nano-black carbon rates on wheat 
crop growth parameters such as days to 
emergence, emergence m7”, plant height, 
spike length, number of spikes m”, 
thousand-grain weight and grain yield. 
Furthermore, we aimed to evaluate the 
suitable level of Nano-black carbon under 
less irrigation to find out the high yielding 
fact that could enhance wheat productivity 
and food security 


2. MATERIALS AND METHODS 
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2.1. Experimental Site 


Wheat performance under _ different 
irrigation levels and various rates of Nano- 
black carbon for grain yield was studied 
under field conditions during the winter 
season of 2018-19 and 2019-20. The soil of 
the experimental site was clay loam with 
less organic matter and slightly acidic. 
Experiments were conducted in Randomized 
Complete Block Design with split plots 
arrangement having three replications. 
Nano-black carbon was allotted to the main 
plots while Irrigation levels were assigned to 
sub-plots. Each experimental unit was 1.5 m 
x 3.3, accommodating 11 rows equally 
spaced at 30cm. Initially, all sowings were 
done at a uniform seeding rate at the rate of 
100 kg ha'. However, the desired 
population was maintained by manual 
thinning. Phosphorous at the rate of 90 kg 
ha! was applied in the form DAP. All the 
phosphorous and % of the nitrogen was 
applied at the time of sowing and Nano- 
black carbon was applied at the time of 
seedbed preparation. 
nitrogen was top-dressed with first irrigation 
and % with second irrigation. The crop was 
sown on a well-prepared seedbed using a 
seed rate of 120 kg ha’!. All other agronomic 
practices were kept normal and uniform for 


The remaining % 


all the treatments 


Following factors were studied during the 
experiment, Nano-black carbon were 
applied to the main plot at the rate of 5 Mg 
ha', 10 Mgha!,15 Mg ha, 20 Mg ha‘'and 
25 Mg ha while irrigation level were treated 
to sub-plots at the rate of 250 mm, 275 mm 
and 300 mm respectively. 
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2.2 ‘Procedure for Recording Data 


Data regarding days to emergence was 
recorded by counting the days taken from the 
date of sowing to the date when 50% 
emergence occur in each plot. Data on 
emergence m7” was recorded in two central 
rows by using the following formula. Data on 
the number of tillers was recorded by 
counting the numbers of tillers in central two 
rows of each plot and was then converted into 
numbers of tillers m?. Data on plant height 
was recorded by measuring randomly selected 
10 plants in each plot from the base of the 
plant to the tip of spikes excluding awns at 
physiological maturity. Data on spike length 
was recorded by using a plastic scale of 
length 30cm.Data regarding grain spike’ was 
recorded by counting wheat grains in 
randomly selected five spikes in each plot and 
was averaged accordingly. Data on grain 
yield was recorded by harvesting the two 
central rows in each plot and was sun-dried, 
weighed, and then converted to kg ha!. Data 
regarding thousand grains weight was 
recorded on the sensitive electronic balance 
after counting a thousand grains for each plot. 


2.3 Statistical Analysis 


The recorded data were statistically 
analyzed according to the analysis of 
variance techniques used for randomized 
complete block design and _ the least 
significant difference (LSD) was used at a 
5% level of significance (P<0.05) upon 
significant F-test through the procedure 
described by Jan et al. (2009). 


3. RESULTS AND DISCUSSION 
3.1 Days to Emergence (m) 
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Data regarding days to emergence (DE) in 
response to different irrigation levels and 
Nano-black carbon showed _ significant 
differences in emergence timing (Table 1). 
Earlier emergence was noted with 5 Mg ha 
followed by 10 Mg ha‘and 15 Mg ha’ 
having at par values while delayed 
emergence was recorded with 25 Mg ha 
followed by 20 Mg ha''respectively (Fig 1). 
The reason for emergence deviation might 
be due to environmental factors which 
influence enzyme activation and rupturing 
of seed coat and ultimately seed emergence. 
The optimum temperature activates most of 
the enzymes in aleoron layer of the seed coat 
which activates the embryo of the seed and 
leads to early emergence. The possibility of 
early emergence in 5 Mg ha’! might be due 
to optimum temperature in field conditions 
and enhanced emergence as compared to 
other Nano-black carbon. The delayed 
emergence might be due to environmental 
stress in which seeds do not react well and 
become dormant due to harsh components of 
the environment (minimum temperatures, 
rainfall, humidity, winds and sunshine, etc) 
and may take more days to emergence. 
Similarly, in the case of wheat Irrigation 
levels, earlier emergence was produced by 
275mm irrigation followed by 250 mm 
irrigation level whereas the delayed 
emergence was noted 300 mm irrigation. 
This might be the genetic character of the 
variety and leading to early or late 
emergence. These results are in connection 
with the findings of Imran et al. (2020) who 
reported that differences in the total roots of 
a young seedling of different wheat 
Irrigation levels were affected by different 
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moisture levels and Irrigation levels. While 
their interactions were non-significant. 
These results are supported by Shah et al. 
(2011) who reported that reduction in 
moisture content reduces the seminal root 
length. The difference in root length is the 
cause of survival of a seedling in stresses 
and its re-growth potentials under abnormal 
situations. The better the root grows with a 
relatively longer length may result inthe 
better establishment of the seedling and 
withstand against abnormal situations e.g. 
drought, high temperatures and salt stress to 
convert the seedling to a healthy plant for 
production (Khan et al., 2014). 


3.3. Emergence (m7) 


Data regarding emergence per m®” of 
different wheat Irrigation levels (Table 2) 
showed significant differences in seedling 
emergence of different wheat Irrigation 
levels as affected by different Nano-black 
carbon. More seedlings were counted with 
seeds sown on 20 Mg ha’ followed by 5 Mg 
ha! and 15 Mg ha’ having at par values 
while less number of emerged seedlings 
were counted with 10Mg ha’! and then by 25 
Mg ha’ respectively (Fig 2). The reason for 
emergence deviation might be due to 
environmental factors which influence 
enzyme activation and rupturing of seed coat 
and ultimately seed emergence. 

The optimal temperature activates the 
majority of the enzymes in the aleoron layer 
of the seed coat, wakes up the embryo of the 
seed, and results in early emergence. When 
compared to other Nano-black carbon rates, 
the prospect of early emergence in 5 Mg ha”! 
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Table 1. Days to emergence as influenced by different Nano-black carbon rates and irrigation 
levels (Data pooled over the both years) 


Nano-black carbon (Mg ha’') Irrigation levels Mean 
250 mm 275 mm 300 mm 
| fhe B 6.333 F335 7.000d 
10 9.000 10.333 9,333 9.556c 
15 9.000 9333 9:333 9272 
20 12.3353 10.000 1333 11.222b 
25 20.667 21.333 24.667 222224 
11.667ab 11.467b 12.400a 
LSD value (P< 0.05) for Irrigation levels = 0.6 
LSD value (P< 0.05) for Nano-black carbon = 0.8 
30 - 
—@— irrigation 250mm 
25 - ~-O- irrigation 275mm 
—-¥— Irrigation 300mm ie 
® 
e 
2. 20:4 
© 
= 
® 
g 
wm 157 
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© 
im) 
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5 10 15 


Nano-Black Carbon rates (Mg ha’) 


Figure 1. Response of days to emergence to different irrigation levels and Nano-black carbon 
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may boost emergence. The emergence may 
be due to environmental stress, in which 
seeds do not react effectively and go latent 
as a result of severe environmental 
components, and it may take additional days 
for the seeds to emerge. Similarly, more 
seedlings were counted with 250mm 
irrigation, followed by 275mm irrigation, 
however minimal emergence (m-2) was 
found with 300mm _ urigation having 
Statistically equal values with 275mm 
irrigation. This might be due to the genetic 
characteristics of the variety, resulting in 
early or late emergence. Different moisture 
levels and Irrigation levels have an impact 
on Irrigation levels. Even though their 
encounters were insignificant. According to 
literature, the root is the most important part 
of a plant and has strived to be the best until 
the availability and search for water to 
support plant development (Ozham and 
Hajibabaei, 2014). The more roots a plant 
has, the healthier it is, the better it can 
withstand drought and wind, and_ the 
seedling has a better chance of survival. Our 
findings are confirmed by the findings of 
Abdoli and Saeidi (2012), who revealed that 
decreasing the moisture level reduces the 
radical weight of different wheat irrigation 
levels. Water stress inhibits the mobilization 
of starchy endosperm in various species 
(Bouaziz and hicks 1990). Increasing early 
growth has the potential to increase soil N 
(Pang et al., 2014) and P (Ryan et al., 2014) 
absorption, hence improving crop nutrient- 
use efficiency and weed competitiveness 
(Coleman et al., 2001). 
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3.4 Number of Tillers (m7) 


Data on the number of tillers per m-2 
revealed substantial changes in the number 
of tillers per m-2 of different irrigation 
levels as impacted by varied Nano-black 
carbon concentrations (Table 3). The 
different Nano-black carbon effect was non- 
significant for the number of tillers m7 
showing statistically at par value for all 
Nano-black carbon. The Irrigation levels 
effect was significant and more tillers were 
counted with 275 mm irrigation followed by 
250 mm irrigation whereas the tinniest tillers 
m” were recorded with 300 mm irrigation 
(Fig 3). Another possible reason could be 
due to environmental factors which 
influence enzyme activation. The optimum 
temperature may enhance plant growth and 
development. The possibility might be due 
to optimum temperature in field conditions 
and enhanced plant lateral growth and 
branches initiations. These results are in 
connection with the findings of Imran et al. 
(2020) who reported that differences in the 
total roots of a young seedling of different 
wheat irrigation levels were affected by 
different moisture levels and _ Irrigation 
levels. While their interactions were not 
significant. 


The possible reason could be increased 
early growth has the potential to increase 
soil N (Pang et al., 2014) and P (Ryan et al., 
2014) absorption, hence improving crop 
nutrient-use efficiency and weed 
competitiveness (Coleman et al., 2001). 
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Table 2. Emergence m” as influenced by different Nano-black carbon rates and irrigation levels 


(Data pooled over the both years). 


Nano-black carbon (Mg ha’') Irrigation levels Mean 

250 mm 275 mm 300 mm 

5 188.33 135.00 139.00 154.1la 

10 77.00 83.67 62.00 74.22b 

15 128.00 105.33 73.67 102.33b 

20 200.67 166.00 158.67 175.1la 

2 104.00 87.33 110.00 100.44b 
139.60a 115.47b 108.67b 


LSD value (P< 0.05) for wheat Irrigation levels = 0.7 


LSD value (P< 0.05) for Nano-black carbon = 0.9 


220 4 —®— Irrigation 250mm 
(e) Irrigation 275mm 
200 + —-¥— Irrigation 300mm 


Emergence (m*) 


Nano-Black Carbon Rates (Mg ha’') 


Figure 2.Response of emergence m”to different irrigation levelsandNano-black carbon 


More robust crops collect more light, 
maximizing crop growth rates and biomass, 
especially with late planting or in settings 
with shorter crop duration (Takahashi and 
Gotoh, 1996; Regan et al, 1997). 
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Appropriate variety selection based on 
seedling performance is thus essential for 
higher yield and crop. growth under 
abnormal conditions, such as drought stress, 
which is most commonly encountered by 
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wheat crop in the country as well as in the 
province where more than 66 percent of the 
area is rain-fed. 


3.5 Plant Height (cm) 


Data regarding plant height (Table 4) 
divulge that the highest plant height was 
recorded with 10 Mg ha! followed by 
statistically similar values of 15 Mg ha’! and 
5 Mg ha. The dwarf plants were produced 
when the wheat Irrigation levels were sown 
on 25 Mg ha’ and 20 Mg ha’ respectively. 
Among the irrigation levels, the highest 
plant height was noted in irrigation level of 
300 mm followed by 250 mm _ having 
statistically similar values. The lowest plant 
height was produced by 275 mm irrigation 
rate respectively. The possibility might be 
due to optimum temperature in field 
conditions and enhanced emergence as 
compared to other Nano-black carbon (Fig 
4). The explanation might be related to 
environmental fluctuations in temperature, 
rainfall, humidity, and so on, which cause 
extra days to develop. Similarly, more 
seedlings were counted with 50 mm 
irrigation level, followed by 275 mm, 
although minimal emergence (m-2) was 
detected with 300 mm irrigation level 
having statistically identical values with 275 
mm. This might be due to the genetic 
characteristics of the variety, resulting in 
early or late emergence. These findings are 
consistent with those of Imran et al. (2020), 
who discovered that variations in the total 
roots of a young seedling of various wheat 
Irrigation levels were influenced by 
moisture and Irrigation levels. Even though 
their encounters were insignificant. Variety 
Hashim-08 has the greatest root number in 
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the group, followed by Pak-2013, DN-84, 
and Pirsabak-2005. While Kpk-2015 has the 
lowest amount of roots reported. It is 
reported that root is the fundamental portion 
of a plant and has striven to the best till the 
availability and search of water to sustain 
plant growth (Ozham and Hajibabaei, 2014). 
When the quantity of roots rises, the crop 
becomes more resilient to drought and wind, 
and the seedling has a better chance of 
surviving even when less water is available. 
More vigorous crops collect more light, 
maximizing crop growth rates and biomass, 
especially when planting late or in settings 
where crop duration is limited (Takahashi 
and Gotoh, 1996; Regan et al. 1997). 
Appropriate variety selection based on 
seedling performance is thus essential to be 
taken into account for higher yield and crop 
growth under abnormal circumstances, such 
as drought stress, which is most commonly 
faced by wheat crops in the country as well 
as in the province where more than 66 
percent of the area is rain-fed. 


3.6 Spike Length (cm) 


Mean values of the data revealed that the 
highest spike length was recorded with 15 
Mg ha’! followed by statistically similar 
values of 5 Mg ha’! and 10 Mg ha’'(Table 
5). Minimum spike length was recorded in 
those plots which were sown on 25 Mg ha’! 
or either 20 Mg ha’. Among the Irrigation 
levels, lengthy spikes were produced with 
300 mm followed by statistically similar 250 
mm and 275 mm irrigation levels (Fig 5). 
The possibility might be due to optimum 
temperature in field conditions and 
enhanced emergence as compared to other 
Nano-black carbon. 
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Table 3. Number of tillers m7as influenced by different Nano-black carbon rates and irrigation 
levels (Data pooled over the both years). 


Nano-black carbon (Mg ha’') Irrigation levels Mean 

250 mm 275 mm 300 mm 

5 200.33 229;33 144.67 191.44 

10 184.00 223.00 140.67 182.56 

15 171.00 233.00 148.33 184.11 

20 191.33 184.00 187.67 187.67 

25 217.33 179.33 169.67 188.78 
192.80ab 209.73a 158.20c 


LSD value (P< 0.05) for wheat Irrigation levels = 10.3 
LSD value (P< 0.05) for Nano-black carbon = 11.2 


—®— Irrigation 250mm 
240 oO. Irrigation 275mm 
Pe —-y¥— Irrigation 300mm 


220 


200 


180 


160 


Number of tillers m° 


140 


120 


5 10 15 20 25 


Nano-Black Carbon rates (Mg ha”) 


Figure 3. Response of number of tillers m™“to different Irrigation levelsandNano-black carbon 
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Table 4. Plant height (cm) as influenced by different Nano-black carbon rates 
levels (Data pooled over the both years). 


Imran. 


and irrigation 


Nano-black carbon (Mg ha’') Irrigation levels Mean 

250 mm 275mm 300 mm 

2 101.83 90.32 94.30 95.48ab 

2 107.21 97.42 112.75 105.79a 

- 90.51 94.95 107.54 97.67ab 

so 92.11 92.11 96.77 93.66b 

= $6.52 80.65 92.72 86.63b 
95.64ab 91.09b 100.82a 


LSD value (P< 0.05) for wheat Irrigation levels = 10.3 
LSD value (P< 0.05) for Nano-black carbon = 11.2 
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Figure 4. Response of plant height to different irrigation levels and Nano-black carbon 
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—e— Irrigation 250mm 
145 ° Irrigation 275mm 
—-¥— Irrigation 300mm 


Spike length (cm) 


Nano-Black Carbon Rates (Mg ha’') 


Figure 5. Response of spike length to 
different Irrigation levels and Nano-black 


carbon. 


The enlarged plant might be attributed to 
temperature, and similarly, more seedlings 
were counted with 250 mm irrigation, 
followed by 275 mm, and_ the least 
emergence (m-2) was reported with 300 mm 
irrigation. These findings are consistent with 
the findings of Imran et al. (2020), who 
discovered variances in the total roots of 
early seedlings of different wheat varieties. 
Different moisture levels and Irrigation 
levels affectedthe irrigation levels. While 
their encounters were insignificant. 
Increasing crop nutrients, vigorous crop 
growth, and biomass have the potential to 
boost crop nutrients, vigorous crop growth, 
and biomass, especially with late planting or 
in conditions where crop duration is shorter. 
Drought stress is the most prevalent problem 
that wheat faces and 66 percent of the area 
are used to support plant development. 


3.6 Thousand Grain Weight (g) 
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Data regarding thousand-grain weight of 
different wheat Irrigation levels as 
influenced by different © Nano-black 
carbon(Table 6). The mean values divulge 
that more thousand-grain weight was 
produced by those plots which were sown 
on 5 Mg ha" followed by 10Mg ha’! whereas 
minimum thousand seed weight was 
recorded with 20 Mg ha’ followed by 15 
Mg ha”! and 25 Mg ha’! respectively (Fig 6). 
The Irrigation levels were non-significant 
for thousand seed weight and noted that all 
the cultivars produced statistically similar 
thousand-grain weight. The fluctuation in 
thousand-grain weight with different Nano- 
black carbon might be due to optimum 
temperature in field conditions and 
enhanced emergence as compared to other 
Nano-black carbon. Plant phenology can be 
improved by using nano-black carbon and 
irrigation. These findings are related to the 
work of Imran et al. (2020), who discovered 
that variances in the total roots of a young 
seedling were impacted by varying moisture 
levels. 


Biochar increases soil N and P to improve 
crop nutrient utilization efficiency and weed 
competitiveness (Coleman et al., 2001). In 
the case of Irrigation levels comparison, the 
highest grain yield was recorded with the 
sowing of the cultivar 250 mm followed by 
300 mm having statistically similar produce. 
More vigorous crops collect more light, 
maximizing crop growth rates and biomass, 
especially when planting late or in settings 
where crop duration is limited (Takahashi 
and Gotoh, 1996; Regan et al., 1997). 
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Table 5. Spike length (cm)as influenced by different Nano-black carbon rates and irrigation 
levels (Data pooled over the both years). 


4 Irrigation levels 
Nano-black carbon (Mg ha”) 250 mm 275mm 500 Gum Mean 
5 10.110 10.233 10.867 10.403ab 
10 9.800 9.710 11.267 10.259ab 
15 10.00 9.843 12.820 10.888a 
20 7.800 9537 10.910 9.416b 
2 (ERE: 7.933 8.267 7.857¢ 
9.129b 9.339b 10.826a 


LSD value (P< 0.05) for wheat Irrigation levels = 10.3 
LSD value (P< 0.05) for Nano-black carbon = 11.2 
Table 6. Thousand grain weight (g) as influenced by different Nano-black carbon rates and 


irrigation levels (Data pooled over the both years). 


Nano-black carbon (Mg ha’') Irrigation levels Mean 

250 mm 275 mm 300 mm 

: 46.667 Da393 46.667 48.889a 

i 46.667 33,393 40.000 40.000b 

» 33.333 30.000 36.667 33.33bc 

a 33.333 33.333 30.000 32.222¢ 

= 36.667 33,333 30.000 33;333de 
39.333a 36.667a 36.667a 


LSD value (P< 0.05) for wheat Irrigation levels = 10.3 
LSD value (P< 0.05) for Nano-black carbon = 11.2 
Table 7. Grain yield (kg ha‘') as influenced by different Nano-black carbon rates and irrigation 


levels (Data pooled over the both years). 


Nano-black carbon (Mg ha’') Irrigation levels Mean 
250 mm 275 mm 300 mm 
5 4404.2 4579.3 4656.2 4546.6 a 
10 3157.6 1777.8 3609.2 3048.2 b 
15 33399 1764.2 1670.1 2258.1 be 
20 2122.6 1631.0 1973.3 1908.9 cd 
2 1710.6 888.7 1333.4 1310.9 d 
3067.0a 2128.2b 2648.4ab 


LSD value (P< 0.05) for wheat Irrigation levels = 10.3 
LSD value (P< 0.05) for Nano-black carbon = 11.2 
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Figure 6. Response of thousand-grain 
weight to different Irrigation levels and 
Nano-black carbon 


Appropriate variety selection based on 
seedling performance is thus essential to be 
taken into account for higher yield and crop 
growth under abnormal circumstances, such 
as drought stress, which is most commonly 
faced by wheat crop in the country as well 
as in the province where more than 66 
percent of the area is rain-fed. Despite the 
fact that their encounters were insignificant. 


3.7. Grain Yield (kg ha“) 


Data regarding grain yield of different 
wheat Irrigation levels as influenced by 
different Nano-black carbon (Table 7). 
Mean data revealed that the highest wheat 
grain yield was produced by those plots 
which were sown on 5 Mg ha’! Followed by 
10 Mg ha" having at par value with 15 Mg 
ha'. The lowest grain yield was produced 
when the plots were sown on 25 Mg ha’! 
followed by 20 Mg ha’ respectively (Fig 7). 
The lowest grain yield was noted with the 
sowing of the cultivar 275 mm. Generally, it 
was concluded that the most productive 
cultivar was 250 mm having the supreme 
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grain yield as compared to other sown 
cultivars. 


Among the Nano-black carbon, either 5 Mg 
ha'for wheat production was noted the 
outmost promising factor to overcome on 
food hunger and inconsistency. 


5000 5 


—*®— Irrigation 250mm 
° Irrigation 275mm 
—-¥— Irrigation 300mm 


4000 5 
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2000 5 
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v 


1000 4 


Nano-Black Carbon Rates (Mg ha”) 


Figure 7. Response of grain yield to 
different Irrigation levelsandNano-black 
carbon 


Imran et al. (2020) stated that plant 
phenology may be enhanced by employing 
Nano-black carbon and watering. They also 
demonstrated that more vigorous crops are 
taken into account for greater yield and crop 
growth in atypical situations. Variety 
Hashim-08 has the greatest root number of 
the bunch, followed by Pak-2013, DN-84, 
and Pirsabak-2005. While Kpk-2015 has the 
fewest roots. According to reports, the root 
is the most important part of a plant and has 
strived to be the best until the availability 
and search for water to support plant 
development (Ozham and Hajibabaei, 2014). 
The bigger the root number, the healthier the 
plant, the better it will be able to withstand 
drought and wind, and the seedling will have 
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a better chance of survival if water is made 
accessible at any stage throughout the crop's 
growth. 


4. CONCLUSION 


According to the findings of this study, a 
correct rate of Nano-black carbon can 
considerably improve the development of 
the root system, which may ultimately boost 
shoot growth and final yield. The current 
investigation found that wheat Irrigation 
levels (250 mm) may germinate properly 
even at low moisture and varying rates of 
Nano-black carbon. Variation in sowing date 
decreases growth and other yield-related 
qualities in a linear fashion. Different 
Irrigation levels responded differently to 
different Nano-black carbon, indicating that 
there is potential for additional development 
in genotype. It is necessary to develop such 
irrigation levels that can be grown in every 
environment and have changeable Nano- 
black carbon capacity to resist diverse 
environmental influences. 
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ABSTRACT 


Mung bean is one of the important Kharif pulses in Pakistan and is grown 
mainly for its edible seeds; therefore, fertilizers management is an 
important factor for improving mungbean growth and yield. A field 
experiment was conducted during the summer of 2013 at Palato Farm of 
the University of Agriculture Peshawar, Amir Muhammad Khan Campus 
Mardan, to determine the effect of phosphorus (P) and Zinc (Zn) on the 
yield and yield component of mungbean. The experiment consisted of four 
levels of P (0, 25, 50, and 75 kg ha‘!) and four levels of Zn (0, 5, 10, and 
15 kg ha‘!). Data associated with the number of leaves and plant height 
illustrated that the higher number of leaves plant"! (8.8) by an average was 
observed when P was applied at the rate of 75 kg ha! followed by 0 kg 
phosphorous (P) ha"! (8.7) and Zn (Zn) application at the rate of 10 kg ha" 
‘produced a maximum number of leaves plant! (9) followed by 15 kg ha 
(8.8) where 0 kg ZN ha'! resulted in (7.7). Similarly, Zn significantly 
affected plant height, while P and interaction between P and Zn levels 
were non-significant. The higher plant height (95.1 cm) was observed 
when P was applied at the rate of 75 kg ha, followed by 50 kg P ha! 
(93.6 cm). Higher plant height (95.8cm) was recorded when ZN was 
applied at the 5 kg ha! followed by 10 kg ha‘(95.1cm). Higher numbers 
of nodules (13.1) were observed with the application of 50 kg P ha’! 
followed by 75 kg P_ ha’! (12.3), while the lowest (10.6) nodules were 
observed in the control plot. P application at the rate of 25 kg ha'! 
produced a higher grain yield than 75 and 50 kg ha-1 and Zn application at 
the rate of 5 kg ha-1 produced a higher grain yield than 10 and 15 kg ha". 
Therefore, a lower rate of P 25 kg ha! and Zn 5 kg ha’! is recommended 
for a higher yield of mungbean in the agro-ecological condition of 
Mardan. 


KEYWORDS: Mungbean, Zinc, Phosphorous, Yield and Yield 
components 
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1. INTRODUCTION 


calcium, 8.5 mg iron, and 49 mg f-carotene 


Mungbean (Vigna radiate L.), also called 
green gram, is an important summer- 
growing pulse crop in Pakistan (Hakim et 
al., 2021; Ali et al., 2019; Ahmad et al., 
2003). It contains 24.5% protein and 59.9% 
carbohydrates. It also contains 75 mg 
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per 100g of split dal (Hakim et al., 2021; 
Shakya et al., 2019). The foliage and stem 
are good sources of fodder for livestock. Its 
seed is more palatable nutritive, digestible, 
and non-flatulent than other pulses (Teferie 
et al., 2020; Tarafder et al., 2020). The 
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unique and common feature of mungbean is 
the root nodules that contain aerobic bacteria 
called rhizobia which fix atmospheric 
nitrogen in the root and thus enhance soil 
fertility (Singh et al, 2021; Ashraf et al., 
2003). It is also a substitute for animal 
protein and forms a balanced diet when used 
with cereals (Detzel et al., 2021). In 
Pakistan, mungbean is cultivated as a minor 
crop and used as food. The area under 
mungbean in Pakistan was around 141 
thousand hectares, with 93 thousand tons in 
2011 (Hakim et al., 2021; GOP, 2012). 


Fertilizers management is one of the 
important factors for improving the growth 
and yield of mung-bean (Ali et al., 2019; 
Iqbal et al., 2021; Ullah et al., 2020). 
Phosphorous is an essential component of 
ADP, ATP, the cell wall, and DNA and 
plays a key role in promoting plants storage 
and structural activities (Aimen et al., 2021). 
P is an important element that significantly 
affects plant growth and _ metabolism 
(Amanullah et al., 2022; Sadiq et al., 2017; 
Bashir et al., 2011) and is a component of 
DNA and RNA, involved in cell division 
(Brady and _ Weil, 2004). Nodule 
establishment and its function are important 
sinks for P, and nodules usually have the 
highest P content in the plant (Sulieman et 
al., 2015). It is supposed that P is effectively 
translocated into grain at high rates since P 
is necessary for producing protein, 
phospholipids and phytin in bean grain 
(Rahman et al., 2008). Poor nodulation and 
poor plant vigour are observed in beans 
grown in P deficient soils (Bindraban et al., 
2020). Among other essential factors, an 
appropriate supply of micronutrients is also 
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required for crops' proper growth and yield. 
The importance of Zn as a micronutrient in 
crop production has increased in recent 
years (Amanullah et al., 2020; Thapa et al., 
2021). Hence Zn is considered the most 
yield-limiting micronutrient (Arunachalam 
et al., 2013). 


The Zn application essentially is being 
employed in the functional and structural 
components of several enzymes (Amanullah 
et al., 2020; Read et al., 2019), such as 
carbonic anhydrase, alcohol dehydrase, 
alkaline phosphatase, phospholipase, 
carboxypeptidase (Read et al., 2019) and 
RNA polymerase (Romheld and Marschner., 
1991). Furthermore, plants that emerged 
from seeds with lower Zn could be 
susceptible to biotic and abiotic stresses 
(Rehman et al., 2018). Zn enriched seeds 
perform better concerning seed germination, 
seedling growth, and yield of crops (Haider 
et al., 2020). In addition, Zn acts as an 
activator of several enzymes in plants and is 
directly involved in the biosynthesis of 
growth substances such as auxin, which 
produces more plant cells (Umair et al., 
2020; Gobarah et al., 2006). 


Furthermore, Zn enhanced photosynthesis 
at the early growth of plants, improved 
nitrogen fixation, grain protein, and yields of 
mungbean plants (Umair et al., 2020). In 
Pakistan, Zn scarcity in the soil is the first 
most widespread problem. In Khyber 
Pukhtunkhwa, the extent of Zn deficient 
soils ranges from 21% to 77%. 42% of 
agricultural fields of Mansehra and Swat 
have a Zn deficiency. On average, 37% of 
fields are deficient in Zn (Ahsin et al., 
2020). 
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Figure 1. Map representing an experimental location in Mardan city, KPK Province, Pakistan, 


using Google maps. 


Therefore, the present experiment was 
designed to study the effect of different 
levels of P and Zn on the yield of mungbean 
and find out the best combinations of P and 
Zn for higher yield and yield components of 
mungbean under the agro-climatic condition 
of Mardan. 


2. MATERIALS AND METHODS 
2.1 Experimental design 


A field experiment was conducted during 
the summer of 2013 at Palato Farm of the 
University of Agriculture Peshawar, Amir 
Muhammad Khan Campus Mardan (Figure 
1), to determine the effect of P and Zn on the 
yield and yield component of mungbean. 
The experiment was conducted in a 
randomized complete block design with 
three replications, and the plot size was 2 m 
x 1.8 m. Mungbean variety “Ramzan” was 
sown in lines having five rows 35cm apart 
on July 2, 2013. The experiment consisted of 
four levels of P (0, 25, 50 and 75 kg ha"!) 
and four levels of Zn (0, 5, 10, and 15 kg ha’ 
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') SSP and ZnSO, will be used as the source 
of P and Zn, respectively and applied as a 
whole during seedbed preparation. All other 
agronomic practices, such as weeding, 
irrigation, plant protection measures, etc., 
were normal and uniform for all the 
experimental units. 


2.2 Data collection and measurement 


The number of plants m® were recorded 
from three central rows by counting the 
number of plants within the metering rod of 
each row and was converted to m7. Data on 
the number of branches plant! was recorded 
by selecting five plants randomly from each 
treatment, and the number of branches were 
counted from base to top and were averaged. 
The number of leaves plant! was recorded 
by selecting five plants randomly from each 
treatment, and the number of leaves was 
counted from base to top and averaged. Data 
on plant height was recorded by randomly 
selecting five plants from each treatment and 
measuring its height from base to tip and 
then averaged to record plant height and 


40 


Journal of soil, plant and Environment 


averaged them. The number of pod plant! 
was recorded by selecting five plants 
randomly from each plot and then were 
picked from it, and the number of pods was 
counted, and then the mean was calculated. 
Data on the number of grain pod! was 
recorded by counting the number of grains 
from five randomly selected pods in each 
plot and then averaged. Numbers of nodules 
plant-1 were recorded by selecting three 
plants randomly from each treatment, and 
the number of nodules were counted on 
roots and averaged. Data on Biological yield 
(kg ha!) was recorded by harvesting three 
central rows in each plot and kept in the 
field for sun drying. It was weighed with the 
help of scale and converted into kg ha’! by 
this formula. 


(Weight of bundles/ No. of rows, *Row 
length * row to row distance) x 10000 


Data on grain yield was recorded by 
harvesting one square meter area from each 
plot and then were threshed, cleaned, dried, 
and weighed. The dried grains were weighed 
with the help of electronic balance and then 
converted to kg ha’!. Data on thousand 
grains weight (g) was recorded by counting 
1000 grains from each plot and then were 
weighed with the help of a sensitive 
electrical balance. The harvest index was 
calculated using the formula. 


Harvest Index (%) = (Grains yield / 
biological yield) x100 


2.3 Statistical analysis 


The data were statistically analyzed using 
analysis of variance techniques appropriate for 
randomized complete block design. Means 
were compared using LSD test at 0.05 level of 
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probability when the F-values were found 
significant (Jan et al., 2009). 


3. RESULTS 
3.1 Phenology and Physiology 


Data recorded on the number of branches 
of plant! are presented in figure 2-A. In 
comparison, P and Interaction between P 
and Zn levels was found non-significant. 
The mean value of the data indicated that a 
higher number of branches plant”! (1.88) was 
recorded when P was applied at the rate of 
75 kg ha! followed by 50 kg P ha’! (1.83) 
where 25 kg P ha’ results number of 
branches plant’ (1.17). The higher number 
of branches plant! (1.93) was recorded 
when Zn was applied at the rate of 10 kg ha” 
' followed by 15 kg ha’!(1.92), .where 0 kg 
Zn ha’! resulted in the lower number of 
branches plant’! (1.58). 


Data recorded on the number of plants m”, 
are presented in figure 2-B. Statistical 
analysis of the data showed that Zn had a 
significant influence on the number of plants 
m” while P and interaction between P and 
Zn levels were found non-significant. The 
mean value of the data indicated that a 
higher number of plants m7? (22.9) was 
recorded when P was applied at the rate of 
75 kg ha'!, followed by a control plot (21.6), 
whereas 25 kg P ha” had the least number of 
plant m®? (20.3). The higher number of 
plants m® (24) was recorded when Zn was 
applied at the rate of 15 kg ha” followed by 
10 kg ha’!(23), where the least number of 
plant m7? (17) were noted in the plot where 
no Zn was applied. 

Data recorded on the number of leaves 
plant’ are shown in figure 2-C. Statistical 
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analysis of the data shows that P and Zn 
significantly influence the number of leaves 
plant’. While Interaction between P and Zn 
levels was found non-significant. The mean 
value of the data indicated that a higher 
number of leaves plant’! (8.8) was observed 
when P was applied at the rate of 75 kg ha’! 
followed by 0 kg phosphorous ha’! (8.7), 
where 50 kg P ha’! resulted in the lower 
number of leaves plant (8). The higher 
number of leaves plant’! (9) was recorded 
when Zn was applied at the rate of 10 kg ha” 
' followed by 15 kg ha"!(8.8), where 0 kg Zn 
ha’! resulted in (7.7) lower number of leaves 
plant’. 

A 


2.5- Zn**. P** _ Zn* 


Branches (Plant) 


ZnO ZnS Zn10 Zn15 


C Zinc (kg ha‘) 


10 | Zn**, P** | Zn* 


b 
L 


Number of leaves (Plant’) 


° 


Zn0 Zn5 2n10 2n15 


Zinc (kg ha‘') 


Number of plants (m) 


Plant height (cm) 
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Data associated with plant height are 
presented in figure 2-D. Statistical analysis 
of the data shows that Zn _ significantly 
affects plant height while P and Interaction 
between P and Zn levels were found non- 
significant. Mean data shows that higher 
plant height (95.1) was observed when P 
was applied at the rate of 75 kg ha’ 
followed by 50 kg P ha’! (93.6), where 0 kg 
P ha’ results from low plant height (88.6). 
Higher plant height (95.8) was recorded 
when Zn was applied at the 5 kg ha’! 
followed by 10 kg ha'(95.1). Where 0 kg Zn 
ha’! resulted in lower plant height (81.9). 


BE Po P2500) Pool p75 


] Zn**, P* , Zn* 


ZnO Zn5 Zn10 Zn15 
D Zinc (kg ha’) 


Fak pk 
95.) Zn", P*,. Zn" 


2) 
o 
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Figure 2. Effect of P and Zn on the number of plants m’, number of branches plant! , number 


of leaves plant’! and plant height (cm) of mungbean. 
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3.2 Grain yield and yield traits (kg ha“) 


Results showed that P and Zn _ levels 
significantly affected the number of pods 
(figure 3A). While Interaction between P 
and Zn levels was found non-significant 
(Figure 3). The mean value of the data 
indicated that higher numbers of pods (10.3) 
were observed when P was applied at the 
rate of 25 kg ha’! followed by 75 kg P ha! 
(10), where the control plot resulted in the 
lower number of pods (9.2). For Zn, higher 
numbers of pods (10.2) were observed when 
Zn was applied at the 5 kg ha! followed by 
10 kg Zn ha™!(10.1). Control plots resulted in 
lower numbers of pods (9). Similarly the 
grains pod” were also affected by Zn and P 
application (Figure 3B). In comparison, P 
and interaction between P and Zn levels 
were found non-significant. Mean data 
shows that a maximum number of grains 
pod! (10.1) was observed when P was 
applied at the rate of 50 kg ha"! followed by 
25 kg P ha! (10). Whereas P application at 
the rate of 50 kg ha! resulted in a lower 
number of grains (9.6). The higher number 
of grain pod! (10.3) was recorded when Zn 
was applied at the rate of 5 kg ha’ followed 
by 10 kg ha‘'!(10), whereas the control plot 
resulted in the lower number of grains (8.7). 


Maximum numbers of nodules (13.1) were 
observed with the application of 50 kg P ha! 
followed by 75 kg Pha’! (12.3), while the 
lowest (10.6) nodules were observed in the 
control plot (Figure 3C). Similarly, higher 
numbers of nodules (12.7) were observed 
with the application of 15 kg Zn ha! and 
followed by the control plot (11.7), while the 
lowest (11.1) nodules were observed on 10 
kg Zn ha’!. Results on grain yield showed 
that that a higher grain yield (826.7 kg ha’') 
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was observed when P was applied at the rate 
of 25 kg ha’, followed by 50 kg P ha’ 
(792.5 kg ha’'), where a lower grain yield 
(737.4 kg ha!) were obtained from the 
control plot. Higher grain yield (835.0 kg ha 
') was recorded when Zn was applied at the 
rate of 5 kg ha’! followed by 15 kg ha’! (800 
kg ha’'). Whereas lower grain yield (719.9 
kg ha’') was recorded at plot had 10 kg Zn 
hal, 


Maximum (4493.4 kg ha’) biological 
yield was observed with the application of 
25 kg P ha"! while the lowest (3543.1 kg ha” 
') biological yield was observed for 75 kg P 
ha’! (Table 1). Similarly, Zn has a maximum 
biological yield (4511.1 kg ha’') with the 
application of 15 kg ha‘!, while the lowest 
(3630.3 kg ha’!) was recorded in the control 
plot. Higher (39.4) thousand-grains weight 
was observed with the application of 25 kg P 
ha’! followed by (38.8g) at the rate of 75 kg 
p ha‘! (Table 2). While the lowest (36.8g) 
thousand grains weight were observed for 
the control plot. Similarly, Zn has a 
significant influence on thousand grains 
weight, having a higher thousand grains 
weight (42.3g) for application of 15 kg Zn 
ha!, while the lowest (34.4g) thousand 
grains weight was recorded in the control 
plot. P and Zn had significant effects, while 
their interaction was non-significant. The 
higher harvest index (21.73%) was observed 
when P was applied at the rate of 75 kg ha! 
followed by 0 kg ha’! (20.54%), where 25 kg 
ha! resulted in a lower harvest index 
(19.30%) (Table 3).. For Zn higher harvest 
index (23.0%) was observed when Zn was 
applied at the rate of 5 kg ha-1, followed by 
a control plot (21.3%). 
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Figure 3. Effects of P and Zn application on pods plant!, grains pod", nodules plant"! and grain 
yield (kg ha!) of mungbean. 


Table 1. Effects of P and Zn on biological yield of mung bean. 


Zn (kg ha‘!) P( kg ha‘!) 
Control 25 50 75 Mean 

Conti 3233.1 4203.9 4301.1 2782.9 3630.3b 
5 3515.4 4456.7 4001.7 3147.7 3780.4b 
10 3885.2 4680.1 3659.4 3476.7 3925.3ab 
15 4328.1 4632.9 4318.2 4765.0 4511.la 

M. 3740.5b 4493.4a 4070.lab 3543.1b 

ean 


Table 2. Effects of P and Zn on thousand-grain weight (gram) of mungbean. 


Zn (kg ha‘!) P( kg ha'!) 
Control 25 50 75 Mean 
Control 29.7 38.0 34.2 35.7 34.4d 
5 37.3 38.0 37.3 37.1 37.4c 
10 38.7 39.3 40.0 39.0 39.3b 
15 41.3 42.3 42.2 43.4 42.3a 
Mean 36.8b 39.4a 38.4a 38.8a 
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Table 3. Effects of P and Zn on harvest index of mung bean. 
Zn (kg ha‘!) P( kg ha‘!) 
control 25 50 75 Mean 
Céntiol 21.2 23.3 20.2 20.5 21.3ab 
5 25.1 16.0 22.7 28.5 23.0a 
10 17.9 18.4 18.1 20.8 18.8b 
15 18.0 19.6 18.5 17.2 18.3b 
M. 20.54a 19.30a 19.89a 21.73a 
ean 


Mean followed by different letters are found significant at 5% level of probabilities. 


Furthermore, the application of 15 kg Zn ha’! 
resulted in a lower harvest index (18.3%). 


4. DISCUSSION 


Micronutrient insufficiency is the primary 
cause of low crop development and 
production in arable soils (Imtiaz et al., 2011). 
Due to intensive agricultural practices, 
unwise use of mineral nutrition, breeding of 
high yielding and advanced varieties, and 
removal of huge quantities of nutrients at 
every crop harvest with lower nutrients 
returns to soils, the degree and extent of 
nutrient deficiency in arable soils has 
recently had serious consequences, resulting 
in lower micronutrients including zinc (Zn) 
in soil (Kanwal et al., 2019). Therefore to 
we conducted a field experiment on different 
Zn rates in combination with P rates to 
determine its effect on mungbean growth, 
yield and yield components. Our results 
showed that P and Zn significantly effected 
number of plants m? , number of branches 
plant! , number of leaves plant! and plant 
height (cm) pods plant! , grains pod", 
nodules plant! and grain yield (kg ha‘') , 
biological yield, thousand grains weight, and 
harvest index. The plots treated with 15 kg 
ha! Zn in combination with 50 kg P ha! 
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Resulted in higher number of plants m”, 


more branches plant! , maximum leaves 
plant’ and maximum plant height (cm), 
higher number of pods plant!, more grains 
pod"',extra nodules plant’ and higher grain 
yield (kg ha’!). The possible explanation for 
these increments might due to the Zn 
promotes nodulation and nitrogen fixation in 
leguminous crops (Masood, et al., 2022; 
Shahrajabian et al., 2022; Gough et al., 
2021). Furthermore Zn is important in the 
formation of auxin, which increases cell 
volume and increases plant height (Wang et 
al., 2016; Oguchi et al., 2004). Cakmak et al. 
(2000) reported that Zn is essential for active 
enzymatic activity, root cell elongation, and 
reducing free radical damage to the cell. 
Another possible explanation for these 
results that when Zn and P applied to soil 
improved soil physical and _ chemical 
properties which consequently improved 
mungbean growth, yield and _ yield 
component (Singh, et al., 2013). A previous 
study documented that compared to control, 
Zn increased branches numbers in plants and 
leaf area of mungbean (Haider et al., 2021). 
Nair studied the genetic diversity of 
mungbean for iron and zinc and discovered a 
large potential for improvement through 
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biofortification, finding 20-40 g Zn 
concentration kg! for dry mungbean seed, 
which was virtually identical in our work. 
Overall the results showed that the addition 
of Zn and P fertilizer to soil can improve 
growth, yield and yield component of 
mungbean. 


5. CONCLUSION 


Our results showed that the growth, yield 
and yield components of mungbean were 
improved in Zn and P fertilizer application. 
On the basis of our results it is concluded 
that P application at the rate of 25 kg ha! 
produced a higher grain yield as compared 
to 75 and 50 kg ha'!. Whereas, Zn 
application at the rate of 5 kg ha’! produced 
a higher grain yield than 10 and 15 kg ha-1 
hence lower rate of 5 kg _ ha-l is 
recommended for higher yield of mungbean 
in agro-ecological condition Mardan. 
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ABSTRACT 


To investigate the response of wheat to different levels of nitrogen (N), 
zinc (Zn) and humic acid (HA), an experiment was conducted at 
Agronomy Research Farm, the University of Agriculture, Peshawar, 
during 2014-15. The experiment was laid out in a randomized complete 
block design having three replications. Three levels of N (80, 120 and 160 
kg ha), Zn (6, 12 and18 kg ha!) and HA (5, 10 and 15 kg ha‘!) were 
used. Results showed that N application at the rate of 160 kg ha’! 
manifested maximum days to physiological maturity (164 days), 
productive tillers m? (248), spikes m7(258), leaf area tiller! (113.6 cm’), 
spike length (10.4 cm), grains spike"! (52), 1000-grain weight (47.5 g), 
biological yield (9260 kg ha‘!), grain yield (3723 kg ha'!) and harvest 
index (40%). Zn treated plots at the rate of 12 kg ha! showed maximum 
days to physiological maturity (162 days), productive tillers m? (241), 
spikes m~ (252), grains spike! (51), 1000-grain weight (45.2 g), biological 
yield (8843 kg ha'!), grain yield (3375 kg ha") and harvest index (39 %). 
Similarly, HA treated plots at the rate of 12 kg ha‘! revealed maximum 
days to physiological maturity (162 days), productive tillers m? (238), 
spikes m7?(249), spike length (9.7 cm), 1000-grain weight (45.00 g), 
biological yield (8649 kg ha"), grain yield (3342 kg ha'!) and harvest 
index (39%). The combined application of N, Zn, and HA had 
significantly affected wheat yield and yield components. It was concluded 
that N at the rate of 160 kg ha'!, Zn 12 kg ha! and HA 10 kg ha 
‘significantly increased yield and yield components of wheat. 


KEYWORDS: Wheat, nitrogen, zinc, humic acid, Tillers, Physiological 
maturity 


This is an open-access article published by the Journal of Soil, Plant and Environment, which permits 
use, distribution and reproduction in any medium, provided the original work is properly cited. 


www.jsape.com 


50 


Journal of soil, plant and Environment 


1. INTRODUCTION 


Wheat (TriticumaestivumL.) is an annual, 
long-day, self-pollinated plant grown in 
winter. It belongs to the family Poaceae 
(Gramineae). It is used as a staple food in 
Pakistan (Ali et al., 2019a). It dominates 
other agronomic crops in production 
(Eeswaran., 2021). Wheat is a rich source of 
carbohydrates and gluten, increasing its 
demand for baking products (Cappelli et al., 
2020). Wheat straw is also used as fodder 
for livestock. It is estimated that 5-10% 
wheat grain is now being consumed as 
poultry and livestock feed (Tricase et al., 
2018). Nitrogen is often the most deficient 
of all the plant nutrients. Nitrogen is the key 
element in achieving consistently high 
yields in cereals (Iqbal et al., 2019; Iqbal et 
al., 2020, Ali et al., 2020; Wu et al., 2021, 
Ali et al., 2022). 


Zinc is an essential micronutrient for 
wheat. Zn plays very important role in plant 
metabolism by influencing the activities of 
hydrogenase and _ carbonic anhydrase, 
stabilization of ribosomal fractions and 
synthesis of cytochrome (Yuvaraj et al., 
2020). Plant enzymes activated by Zn are 
involved in carbohydrate metabolism, 
maintenance of the integrity of cellular 
membranes, protein synthesis, and 
regulation of auxin synthesis and pollen 
formation (Gondalet al., 2021). Zn is a 
micronutrient that is required for plant 
growth relatively in a smaller amount. Zn 
plays a vital role in the physiological 
process of wheat plants such as cell 
elongation, cell 
translocation, 


maturation, sugar 
meristematic tissues 
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development and _ protein 
(Mumivand et al., 2021). 


synthesis 


Nitrogen and Zn together play a 
significant role in crop production. The 
importance of nutrients (micro and macro) 
for the normal growth of crop plants is 
universally recognized. Zn is an essential 
element present in plant enzymatic systems 
(Yuvaraj et al., 2020). Zn has vast numbers 
of functions in plant metabolism and 
consequently, Zn deficiency has many 
effects on plant growth (Umair et al., 2020). 
Zn deficiency is a worldwide nutritional 
constraint for crop production in many types 
of soil worldwide, particularly in cereals 
growing on calcareous soil (Rehman et al., 
2020; Amanullah, 2020). In Pakistan, Zn 
deficiency is a general micronutrient 
disorder on calcareous soil and is considered 
the third most common deficient nutrient 
after N and phosphorous (Amanullah, 2020). 


Humic acid is complex organic molecules 
that are formed by the breakdown of organic 
matter. HA influences soil fertility through 
its effect on the soil's water-holding capacity 
(Ali et al., 2019b; Sutton and Sposito, 2005). 
HA constitutes a stable fraction of carbon 
that improves soil characteristics such as 
improved water holding capacity, pH 
buffering, and thermal insulation (Izhar et 
al., 2020). It also increases N use efficiency 
and stimulates shoot and root growth (Leite 
et al., 2020; Lodhiet al., 2013). HAs 
contributed to soil stability and soil fertility, 
leading to exceptional plant growth and 
micronutrient uptake. HA is an organically 
charged bio-stimulant that significantly 
affects plant growth and development and 
increases crop yield. It has been extensively 
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investigated that HA _ improves - soils’ 
physical, chemical, and biological properties 
(Zaremanesh et al., 2020; Karim et al., 2020; 
Skowronska et al., 2020). 


Nitrogen, Zinc and HumicAcid are 
essential nutrient which play a pivotal role 
in increasing fertility of the soil which in 
turn increases the yield of the crop. Keeping 
in view the above facts and figures, the 
present investigation was _ therefore 
conducted to determine the best level of N, 
Zn and HA for obtaining maximum yield 
and yield components of wheat. 


2. MATERIALS AND METHODS 


2.1 Field location and_ experiment 
management 

The experiment on the response of wheat 
to different levels of N, Zn and HA on the 
yield and yield components of wheat was 
conducted at Agronomy Research Farm, The 
University of Agriculture, Peshawar, 
Pakistan, during winter 2014-15. The 
research comprised of three factors (1) N N 
kg ha! (NI 80 kg ha!, N2 120 kg ha™ and 
N3 160 kg ha") (ii) Zn kg ha‘! (Zn: 6 kg ha", 
Zn 12 kg ha’, Zn3 18 kg ha’) and (iii) HA 
kg ha! (HA; 5 kg ha!, HA? 10 kg ha’, HA 
15 kg ha! )and one control plot for each 
treatment where no N, Zn and HA applied. 
All three factors were used in the following 
combinations. T;: (Control), T2: NiZniHAj, 
T3: Ni Zn, HA2, Ta: Ni Zn HAs, Ts: 
NiZn2HAy, Te: N1Zn2HAo, T?: NiZn2HAs3, Ts: 
NiZn3HAy, To: NiZn3HAz,T 10: Ni Zn3HAs, 
Tu: NoZn HA, Ti2: NoZn HA, T13: 
No2Zn1 HAs, T14: NoZmH Ai, Tis: NoZnm2HAd, 
Tie: N2Zn2HAsz, T17: No2Zn3HAy, Tis: 
No2Zn3H Ao, Tho: N2Zn3HAs, T20: N3Zm HA, 
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Tai: N3ZnHAg, T22: N3Zni HAs, T23: 
N3Zn2H Ay, T24: N3ZmHAd, Tos: N3ZmHAs3, 
T2.6: N3Zn3HAy, T27: N3Zn3HAg, T2s: 
N3Zn3HA3. The experiment was laid out in a 
randomized complete block design having 
three replications. Wheat variety ‘Atta 
Habib-2010’ was sown on 20" November 
2014 at the rate of 120 kg ha" in a plot size 
of 3 x 1.8 m having 6 rows 30 cm apart. At 
the time of sowing full dose of Zn, HA, and 
half dose of N was used while the rest of 
50% N was applied at second irrigation. A 
basal dose of phosphorous (P20s) 100 kg ha’! 
was applied to each plot. Urea, diammonium 
phosphate (DAP), Zn sulphate and HA was 
used as a source for N, P, Zn and HA 
respectively. Irrigation was done according 
to the need of crop. Weedicides was used for 
control of common weeds of wheat. The 
crop was harvested on May 15, 2015. 
Excluding the treatments, all the further 
agronomic practices were done normal and 
uniform. 


2.2 Data collection and analysis 


Data concerning days to emergence was 
recorded by counting the number of days 
from sowing to 80% of the seedlings 
emerged in each plot. Emergence m” data 
was recorded by measuring the number of 
plants emerged in one meter row length at 
three randomly selected rows in each plot 
and was converted to emergence m” using 
the following formula (1); 


: Total number of seedling emerged 
Emergence m*= ————_———— x Im 
Row - rowdistancex row lengthx No.of rows 
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Data on days to anthesis was recorded 
when 80% of plants extruded the anthers in 
each treatment. Days were counted from the 
date of sowing till the date of anthers 
extrusion in each plot. Days to physiological 
maturity were recorded by measuring the 
days from date of sowing to the date when 
80% plants get physiologically mature 
(Anderson et al., 1985). Leaf area tillers”! 
was calculated at the anthesis stage from the 
leaves of ten tillers in four central rows 
selected randomly at each plot. Leaf area 
tiller! was calculated by multiplying the 
average leaf length (cm), leaf width (cm), 
number of leaves, and correction factor then 
dividing by ten tillers. 


Leaf area tiller'= No. leaves x Avg Leaf 
width (cm) xAvg leaf length x CF/ No. tiller 


Correction factor = 0.75 (Abbas et al., 2014) 


For spike length data, the length of ten 
randomly selected spikes was measured 
from the basal joint of the spike till the top 
of the spike and then averaged. Spikes m? 
data was recorded by counting the number 
of spikes in central four rows of each plot 
and then were converted into spikes m”. The 
data on the number of productive tillers m? 
was recorded by counting the tillers having a 
spike in place one meter selected at three 
time randomly in each plot and were 
converted into productive tillers m*. Ten 
randomly selected spikes for grains spike-1 
of wheat were estimated and then averaged. 
Similarly, a thousand grains were randomly 
taken from each plot and the weight was 
calculated with the help of electronic 
balance. Whereas, the biological yield was 
recorded by harvesting four central rows in 
each plot; tied into bundles, sun-dried, 
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weighed and then biological yield was 
calculated in kg ha” by using formula ; 


Hidlogjcal yi Gig ha”) = Biological yield in four central rows 


The four central rows was harvested, dried, 
threshed, cleaned, weighed and_ then 
converted into kg ha!. For Harvest index 
was calculated by using the following 
formula; Harvest index (%) = 
Economic yield 


x 100 
Biological yield 


2.3 Statistical Analysis 


Analysis of variance procedure was 
followed for the statistical analysis of 
recorded data according to the design used. 
Means compared using least significant 
differences (LSD) test at P < 0.05 upon 
significant F-test (Steel, 1997). 


3. RESULTS 
3.1 Phenology and Growth traits 


Analysis of the data on days to 
emergence and emergence m® showed that 
N, Zn, HA, control vs. rest and all the 
interaction between the nutrients were found 
non-significant (Figure. 1). Days to anthesis 
were significantly affected by N and control 
vs. rest had (Figure. 2). While Zn, HA and 
all the interaction were non-significant for 
days to anthesis. The treated plots took more 
days to anthesis (134) as compared to 
control (133). N's mean data revealed that 
minimum days to anthesis (133) were 
recorded at 80 kg N ha’. Days to anthesis 
increased with each increment of N and 
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maximum days to anthesis (134) was 
observed at 160 kg N ha. Similarly data on 
days to physiological maturity revealed that 
N, Zn, HA, control vs. rest and N x Zn 
interaction significantly affected days to 
physiological maturity, while all other 
interactions were non-significant (Figure 
2B). The treated plots showed more days to 
physiological maturity (163) as compared to 
control (156). N's mean data exhibited that 
minimum days to physiological maturity 
(161) were noted at 80 kg N ha’. Days to 
physiological maturity increased with each 
increment of N and maximum days to 
physiological maturity (165) was noted at 
160 kg N_ ha'. Minimum days to 
physiological maturity (162) were observed 
at 6 kg Zn ha’. Days to physiological 
maturity increased with each increment of 
Zn and full days to physiological maturity 
(163) was observed at 18 kg Zn ha. Low 
level of HA (5 kg ha!) resulted in 
minimum days to physiological maturity 
(162). Days to maturity increased with each 
increment of HA The maximum days to 
physiological maturity (163) were revealed 
at 15 kg HA ha’. Interaction among N x Zn 
showed that days to physiological maturity 
increased with each increment of N up to 
160 kg N ha! with 6 kg Zn ha. Similar 
increase in days to physiological maturity 
was observed with each increment of N with 
12 and 18 kg Zn ha’. Data regardingleaf 
area tiller! revealed that N, Zn, HA, control 
vs. rest, N x Zn and N x Zn x HA 
interaction had significantly affected leaf 
area tiller'(Figure 3A), while all other 
interactions were non-significant for leaf 
area tiller’. The treated plots showed 
maximum leaf area tiller! (111.01) as 
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compared to control (99.00). Mean data for 
N showed that minimum leaf area tiller! 
(108.84) was observed at 80 kg N ha’!. Leaf 
area tiller! increased with each increment of 
N and higher leaf area tiller'(113.66) was 
observed at 160 kg N ha’. Minimum leaf 
area tiller’ (108.88) was noted at 6 kg Zn 
ha"'. Leaf area tiller! increased with each 
increment of Zn and maximum leaf area 
tiller! (112.63) was recorded at 18 kg Zn ha” 
' Lower level of HA (5 kg ha’) resulted in 
lowest leaf area tiller! (107.88). Leaf area 
tiller! increased with HA increment and 
higher leaf area tiller! (112.44) was 
recorded at 15 kg HA ha’. Interaction 
between N x Zn revealed that leaf area tiller” 
' increased with each increment of N up to 
160 kg N ha! with 6 kg Zn ha. Similar 
increase in leaf area tiller’ was noted with 
each increment of N with 12 and 18 kg Zn 
ha’!. In the case of Plant height lower level 
of HA (5 kg ha'') produced in smaller plant 
height (88.3), plant height increased with the 
increment of HA and taller plant height 
(92.6) was obtained at 15 kg HA ha’! (Figure 
3B). Interaction between N x Zn showed 
that plant height improved with each 
increment of N up to 160 kg N ha” with 6 
kg Zn ha. Alike increase in plant height 
was noted with each increment of N with 12 
and 18 kg Zn ha’. Data on spike length of 
the data pertaining that N, Zn, HA, control 
vs. rest and N x Zn x HA interaction had 
significantly affected spike length (Figure 
4A). While all other interactions were non- 
significant spike length. The treated plots 
resulted in greater spike length (9.7) than 
control (8.1). Mean data for N showed that 
smaller spike length (9.2) was recorded at 80 
kg N ha’. Spike length enhanced with each 
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Figure 1. Changes in Days to emergence and emergence (m7) of wheat to different N, Zn and HA 
applications. T1 to T28 indicates different treatments, for detail see material and methods. 
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Figure 2. Changes in Days to anthesis and Daysof wheat to different N, Zn and HA applications. T1 to 
T28 indicates different treatments, for detail see material and methods. 
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Figure 3. Changes in leaf area tiller’ and plant height of wheat to different N, Zn and HA applications.T1 
to T28 indicates different treatments, for detail see material and methods. 
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Figure 4. Changes in leaf area tiller-1 and wheat plant height to different N, Zn and HA applications.T1 
to T28 indicates different treatments, for detail see material and methods. 
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increment of N and greater spike length 
(11) was observed at 160 kg N ha’. 
Minimum spike length (9.2) was recorded at 
6 kg Zn hal. Spike length enhanced with 
each increment of Zn and maximum spike 
length (10) was obtained at 18 kg Zn ha. 
The lower level of HA (5 kg ha’) resulted 
minimum spike length (9.1). Spike length 
increased with each increment of HA and 
maximum spike length (9.9) was recorded at 
15 kg HA ha’, 


The application of N, Zn, HA, control vs. 
rest and Zn x HA interaction significantly 
affected spikes m” of wheat (Figure 4B), 
while all other interactions were non- 
treated plots showed 
maximum spikes m” (249) as compared to 
control (209). N's mean data revealed that 
minimum spikes m-2 (239) were noted at 80 
kg N ha’. Spikes m? improved with each 


significant. The 


increment of N and maximum spikes m~ 
(258) was recorded at 160 kg N ha’. 
Minimum spikes m°* (238) was recorded at 6 
kg Zn ha!. Spikes m? improved with each 
addition of Zn and maximum spikes m7” 
(253) was obtained at 18 kg Zn ha‘!. Lower 
level of HA (5 kg ha”) produced minimum 
spikes m” (238). Spikes m® increased with 
each increment of HA and more spikes m7” 
(253) was noted at 15 kg HA ha’. Interaction 
between Zn x HA showed that spikes m7 
enhanced with each increment of Zn up to 18 
kg N ha! with 5 kg HA ha’. Alike increase 
in spikes m* was noted with each increment 
of Zn with 10 and 15 kg HA ha’. 


3.2 Yield Component and Yield Traits 
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The addition of N, Zn and HA 
significantly affected productive tillers m™ 
(Table 9), whereas the interaction among Zn, 
N and HA were found non-significant. The 
treated plots exhibited maximum productive 
tillers m® (239) as compared to control (187). 
Mean data for N exhibited that minimum 
productive tiller m* (230) was recorded at 80 
kg N ha’. Productive tillers m? increased with 
each increment of N and maximum productive 
tillers m? (248) was noted at 160 kg N ha’. 
Minimum productive tillers m* (232) was 
obtained at 6 kg Zn ha’. Productive tillers m 
improved with each increment of Zn and 
higher productive tillers m7 (243) was 
recorded at 18 kg Zn ha’'. Lower level of HA 
(5 kg ha'') resulted lower productive tillers m™ 
(232). Productive tillers m? increased with 
each increment of HA and maximum 
productive tillers m? (242) was noted at 15 kg 
HA ha’. Interaction between Zn x HA 
revealed that productive tillers m* increased 
with each increment of Zn up to 18 kg N ha! 
with 5 kg HA ha’. A similar increase in 
productive tillers m? was recorded with each 
increment of Zn with 10 and 15 kg HA ha’. 


Analysis of the data showed that N, Zn, 
HA, control vs. rest and N x HA interaction 
had significantly affected number of grains 
spike! (Table 11). At the same time, all 
other interaction were non-significant. The 
treated plots revealed more grains spike’! 
(51) as compared to control (36). Mean data 
for N manifested that minimum grains spike” 
' (48) was noted at 80 kg N ha’. Grains 
spike enhanced with each addition of N 
and maximum grains spike! (52) was 
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Figure 5. Changes in productive tillers mand grains spike"leaf area tiller-1 and plant height of wheat to 
different N, Zn and HA applications.T1 to T28 indicates different treatments, for detail see 


material and methods. 


recorded at 160 kg N ha’. Lesser grains 
spike (47) was observed at 6 kg Zn ha’. 


Grains spike’ improved with each 
increment of Zn and more grains spike” (51) 
was obtained at 18 kg Zn ha’'. Lower level 
of HA (5 kg ha’) produced in lesser grains 
spike’ (47). Grains spike! increased with 
each increment of HA and more grains 
spike’ (51) was obtained at 15 kg HA ha’. 
Interaction between N x HA showed that 
grains spike’! improved with each increment 
of N up to 160 kg N ha! with 5 kg Zn ha’. 
Alike increase in grains spike! was recorded 
with each increment of N with 10 and 15 kg 
Zn ha’. Data on 1000 grains weight 
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revealed that N, Zn, HA, control vs. rest, N 
x Zn, Zn x HA and N x Zn x HA interaction 
had significantly affected 1000-grains 
weight (Table 1). While other interaction 
were non-significant. The treated plots 
showed more 1000-grains weight (45) as 
compared to control (36). Mean data for N 
showed that less 1000-grains weight (41.3) 
was obtained 1000-grains weight improved 
with each increment of N and higher 1000- 
grains weight (47.5) was noted at 160 kg N 
ha!. Minimum 1000-grains weight (42.4) 
was recorded at 6 kg Zn ha’. 1000-grains 
weight increased with each increment of Zn 
and maximum 1000-grains weight (45.7) 
was recorded at 18 kg Zn ha’'. Lower level 
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of HA (5 kg ha’) resulted in minimum 
1000-grains weight (43.0). 1000-grains 
weight increased with the increment of HA 
and 1000-grains weight maximum 1000- 
grains weight (45.5) was noted at 15 kg HA 
ha’. Interaction between N x Zn manifested 
that 1000-grains weight increased with each 
increment of N up to 160 kg N ha” with 6 
kg Zn ha'. Similar improved 1000-grains 
weightwas noted with each increment of N 
with 12 and 18 kg Zn ha’. Interaction 
among Zn x HA showed that 1000-grains 
weight increased with each increment of Zn 
up to 18 kg N ha’ with 5 kg HA ha". Alike 
increase in 1000-grains weight was observed 
with each increment of Zn with 10 and 15 
kg HA ha’. 


Analysis of the data revealed that N, Zn, 
HA, control vs. rest, N x Zn and N x Zn x 
HA interaction had significantly affected 
biological yield (Table 2). While all other 
interaction were non-significant. The treated 
plots showed maximum biological yield 
(8805 kg ha!) as compared to control (6127 
kg ha’). Mean data for N revealed that 
lower biological yield (7450 kg ha) was 
noted at 80 kg N ha’. Biological yield 
increased with each increment of N and 
maximum biological yield (9260 kg ha’) 
was observed at 160 kg N ha!. Minimum 
biological yield (7843 kg ha‘') was noted at 
6 kg Zn ha". Biological yield increased with 
each increment of Zn and maximum 
biological yield (8973 kg ha‘') was noted at 
18 kg Zn ha’!. The lower level of HA (5 kg 
ha!) resulted in lowest biological yield 
(7997 kg ha’'). Biological yield increased 
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with increment of HA and higher biological 
yield (8900 kg ha!) was observed at 15 kg 
HA ha’. Interaction between N x Zn 
revealed that biological yield increased with 
each increment of N up to 160 kg N ha! 
with 6 kg Zn ha. Similar increase in 
biological yield was noted with each 
increment of N with 12 and 18 kg Zn ha’. 


The data analysis showed that N, Zn, HA, 
control vs. rest, N x Zn and N x Zn x HA 
interaction had significantly affected grain 
yield (Table 3). While all other interaction 
were non-significant. The treated plots 
showed maximum grain yield (3420 kg ha") 
as compared to control (2360 kg ha”). N's 
mean data revealed that lower grain yield 
(2862 kg ha’') was observed at 80 kg N ha’. 
Grain yield increased with each increment of 
N and maximum grain yield (3723 kg ha’!) 
was observed at 160 kg N ha!. Minimum 
grain yield (3080 kg _ ha’!) was noted at 6 
kg Zn ha’. Grain yield increased with each 
increment of Zn and maximum grain yield 
(3544 kg ha’') was noted at 18 kg Zn ha’. 
Lower level of HA (5 kg ha”) resulted in 
lowest grain yield (3187 kg ha”). Grain 
yield increased with the increment of HA 
and higher grain yield (3493 kg ha’') was 
observed at 15 kg HA ha’. Interaction 
between N x Zn revealed that grain yield 
increased with each increment of N up to 
160 kg N ha! with 6 kg Zn ha. Similar 
increase in grain yield was noted with each 
increment of N with 12 and 18 kg Zn ha 
' Analysis of the data showed that N, Zn, 
HA, control vs. rest, N x Zn and N x Zn x 
HA interaction had significantly affected 
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Table 1. 1000-grain weight (g) of wheat as affected by nitrogen (N), Zinc(Zn) andhumic 


acid(HA). 
N (kg ha’) 
Zn HA Mean 
ats?) gone 80 120 160 
N x Zn x HA Zn x HA 

6 =) 39.9 40.3 42.3 41.5 
10 41.6 42.8 42.8 41.8 

15 42.2 41.5 42.5 42.1 

12 ) 44.5 45.9 46.4 45.6 
10 46.3 45.4 45.1 45.6 

15 46.9 44.6 46.9 46.1 

18 5 44.3 47.1 48.7 46.7 
10 45.3 47.8 49.1 47.6 

15 47.3 48.6 49.6 48.3 

NxHA 
S) 40.0 43.3 45.6 43.0 b 
10 41.8 45.6 47.6 45.0 a 
15 42.1 46.1 48.3 45.5a 
Nx Zn 
6 39.8 43.0 44.4 42.4b 
12 42.2 45.5 47.9 45.2a 
18 41.9 46.0 49.1 45.7a 
41.3¢ 45.8 b A7.5a 


Planned Mean Comparison 


Control 36 
Rest 45 


Note: Means of same category followed by different letters are significantly different at (P< 0.05) using 
LSD test. 


www.jsape.com 60 


Journal of soil, plant and Environment Iqbal et al. 


Table 2. Biological yield (kg ha’) of wheat as affected by nitrogen (N), Zinc(Zn) andhumic 


acid(HA). 
N (kg ha”) 
Zn HA Mean 
(ke ha!) (ke ha!) 80 120 160 
N x Zn x HA Zn x HA 

6 5 7610 7724 7738 7691 

10 7634 7631 8150 7805 

15 7492 7818 8200 7837 

12 5 8432 8850 9520 8934 

10 8915 9282 9394 9197 

15 9011 9152 9671 9278 

18 5 9099 9604 9739 9481 

10 9259 9730 9866 9618 

15 9562 9864 9786 9737 

Nx HA 
> 7230 7880 8880 7997 b 
10 7636 8660 9650 8649 a 
15 7685 9278 9737 8900 a 
Nx Zn 
6 7120 8160 8250 7843 b 
12 7620 9178 9733 8843 a 
18 7609 9528 9797 8978 a 
7450 c 8955 b 9260 a 


Planned Mean Comparison 


Control 6127 
Rest 8805 


Means of same category followed by different letters are significantly different at (P< 0.05) using LSD 
test. 
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harvest index (Table 4). While all other 
interaction were non-significant. The treated 
plots showed a maximum harvest index (39 
%) compared to control (37 %). Mean data 
for N revealed that lower harvest index (38 
%) was observed at 80 kg N ha”. 


Harvest index increased with each 
increment of N and maximum harvest index 
(40 %) was observed at 160 kg N ha’. 
Minimum harvest index (38 %) was noted at 
6 kg Zn ha’. Harvest index increased with 
each increment of Zn and maximum harvest 
index (39 %) was noted at 18 kg Zn ha”. 
Lower level of HA (5 kg ha") resulted in 
lowest harvest index (38 %).Harvest index 
increased with increment of HA and higher 
grain yield (39 %) was observed at 15 kg 
HA ha‘. Interaction between N x Zn 
revealed that harvest index increased with 
each increment of N up to 160 kg N ha! 
with 6 kg Zn ha’. Similar increase in 
harvest index was noted with each increment 
of N with 12 and 18 kg Zn ha’. 


4. DISCUSSION 


Our results showed that days to emergence 
and emergence m7 of wheat to N, Zn, HA, 
control vs. rest and all the interactions 
among the nutrients were non-significantly 
affected. The possible reason may be that 
the consumption of its own stored food in 
the seed as an enough source for seedling to 
emerge. Our results are similar to (Waraich 
et al., 2007) who reported that N and HA did 
not affect emergence. Similar results are 
also reported by Leghari et al. (2016). It may 
be because fertilizer response is not so quick 
and observed immediately after two weeks. 
Similarly, Leghari et al. (2016) reported that 


www.jsape.com 


Iqbal et al. 


emergence m~ was not affected by N. The 
reason may be that seed used its own 
endosperm food for germination and plants 
did not use nutrients from outside source 
(Shah et al., 2012). The results also conform 
to Shah et al. (2009) who stated that 
fertilizer fertilization had no influence on the 
emergence m”. 


N, Zn, HA, control vs. rest, and N x Zn, 
wheat's days to anthesis and maturity were 
significantly affected, while the early 
anthesis was recorded inin control plots. The 
delayed anthesis is due to the nutrients 
available in an adequate amount, which 
enlarged the growing time period of cereal 
crop (Arif et al, 2006). N_ fertilization 
increased vegetation, leaf area and light use 
efficiency (Zeidan et al., 2010). Similarly, N 
application delayed days to tasseling, silking 
and maturity in maize (Amanullah et al., 
2008&Arif et al, 2006). However, the 
application of Zn and HA have no impact on 
days to anthesis (Nawab et al., 2011). 
Physiological maturity delayed with 
increased levels of N, Alike results were 
stated by Sivasankar et al. (1993) who stated 
that N is very imperative for crop 
development and growth. Fertilization of N 
promotes lived green foliage duration 
(Frederick &Camberato, 1995). Increase in 
N levels significantly delayed days to 
maturity Ayoub et al. (1994). The delayed in 
days to maturity may be due to enough 
nutrient accessibility, which finally 
enhanced the growing time period of cereal 
crops Ullah et al. (2021). Alike results were 
reported by Zeidan et al. (2010) who 
reported that N _ fertilizer improves 
vegetative growth, increased light use 
efficiency and grain filling period. Higher 
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levels of Zn and HA increased the time of stated that Zn and HA application with 
maturity of wheat crop. Comparable results higher rates late maturity of the crop which 
were reported by Rajput et al. (2004) who due to easily availability of 


Table 3. Grain yield (kg ha") of wheat as affected by nitrogen (N), Zinc(Zn) and humic acid(HA). 


Zn HA ae Mean 
(kg ha‘!) (kg ha‘!) Nd NA ae 
Nx Zn x HA Zn x HA 

6 5 2933 3060 3143 3045 

10 2886 2930 2970 2929 

15 2792 3153 3233 3059 

12 5 3167 3503 3624 3431 

10 3415 3485 3498 3466 

15 3511 3355 STIS 3547 

18 5 3337 3797 3860 3665 

10 3497 3923 3969 3796 

15 3800 4057 4049 3969 

NxHA 
5 2710 3340 3510 3187 b 
10 2890 3425 3710 3342 a 
15 2956 3547 3975 3493 a 
Nx Zn 
6 2732 3214 5295 3080 b 
12 2934 3395 3195 3375 a 
18 2919 3632 4080 3544 a 
2862 c 3414b 3723 a 
Planned Mean Comparison 
Control 2360 
Rest 3420 


Means of same category followed by different letters are significantly different at (P< 0.05) using LSD 
test. 


www.jsape.com 63 


Journal of soil, plant and Environment Iqbal et al. 
Table 4. Harvest index (% ) of wheat as affected by N, Zn and HA. 
N (kg ha’) 
Zn HA Mean 
(kg ha) (kg ha!) 80 120 160 
N x Zn x HA Zn x HA 
6 5 38 a7 39 38 
10 38 38 39 38 
15 ot 39 39 38 
12 5 38 38 38 38 
10 38 38 | 38 
15 39 a) 39 38 
18 2) 37 40 4] 39 
10 38 40 42 40 
15 40 4] 42 4] 
NxHA 
5 38 38 39 38 b 
10 38 38 40 39a 
15 38 38 4] 39a 
N x Zn 
6 38 38 38 38 b 
12 39 38 40 39a 
18 38 38 42 39a 
38 c¢ 38 b 40a 
Planned Mean Comparison 
Control 62) 
Rest 39 


Means of same category followed by different letters are significantly different at (P< 0.05) using LSD 


test. 


essential nutrients thus enhanced vegetative 
growth of crop which delayed maturity. Leaf 
area tiller! and plant height had 
significantly affected by N, Zn, HA, control 
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vs. rest, Nx Zn and Nx Zn x HA 
interactions. Whereas, all further interaction 
were non-significant. Leaf area tiller! 
enhanced with each accumulation of N and 
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maximum leaf area tiller’! was observed at 
160 kg N ha!. While minimum leaf area 
tiller! was noted in control plots. Our results 
conform with Vogeler et al. (2020) who 
stated that N fertilizer increased leaf area of 
wheat. This may be because Zn and HA 
fertilization increased leaf area due to 
increased vegetative growth. Our results are 
in contract with Arif et al. (2006) and Iqbal 
et al. (2020) who stated that significant 
increase was found in the height of crop 
with fertilization of N. Increased plant 
height due to enough amount of N may be 
credited to additional vegetative growth and 
improvement as a result in intermodal 
extension. The soil application of Zn 
increased plant height linearly. Plant height 
increased 18.7 % over control by 
fertilization of Zn (Ali et al., 2013). 


Additionally, spikes m°, tillers m7, spike 
length, grains spike’ and 1000-grains 
weight revealed that N, Zn, HA, control vs. 
rest and Zn x HA interaction had 
significantly affected. Whereas all other 
interactions were non-significant. The 
nutrients treated plots manifested more 
spikes m? as compared to control. Spikes m™ 
* increased with each increment of N and 
maximum spikes m* was noted at 160 kg N 
ha’. The reason may be that the availability 
of nitrogen in enough amounts facilitates the 
plant's tillering ability, which is directly 
related to spikes m*. Our results are in line 
with Jan et al. (2000) who find that spikes 
population and number of grains spike"! had 
a significantly increased with increasing 
levels of N fertilizer. Alike results were 
published by Shahab et al. (2016). Spikes m 
increased with each addition of Zn and 
maximum spike m* was observed at 18 kg 
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Zn ha’'. Our results are in conformity with 
Keram et al. (2014) who postulated that Zn 
at the rate of 20 kg ha! along with a 
recommended dose of NPK improved yield 
and yield-related components of wheat. 
Spikes m? increased with each incensement 
of HA and maximum spikes m” was noted 
at 15 kg HA ha’. The reason might be that 
HA stimulates microbial movement in soil 
to improve the physical structure and 
biological environment of soil, which 
enhances plant growth and development 
(Zancani et al., 2009). The possible reason 
increase in spike length might be due to that 
N fertilization increased the vegetative and 
reproductive growth of plant (Shahab et al., 
2016). Our results are in agreement with 
Leghari et al. (2016). Zn fertilization also 
increased spike length linearly with each 
increment of Zn and maximum spike length 
was noted at 18 kg ha! Zn. Alloway et al. 
(2004) stated that the combine resulting 
from the effects of N application helps 
enhance plant growth, change the pH of the 
root environment, and found significant 
interactions increasing dose of Zn and N. Zn 
plays a vital role in protein degradation and 
protein-protein interactions. This way, it 
takes part in plant metabolism, growth, and 
development (Marschner et al., 1995). Spike 
length increased with each increment of HA 
and maximum spike length was noted at 15 
kg ha' HA. Grains spike! might be 
improved because Zn's application makes 
readily available nutrients uptake from the 
soil to plant and enhances their effectiveness 
greatly in the grains creation time and 
resulted maximum grains spike! (Zeidan et 
al., 2010). HA application promotes plant 
growth by providing unavailable nutrients 
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and buffering pH of soil (Julie and Bugbee, 
2006). Improvement in 1000-grains weight 
may be due to that plant protein synthesis, 
production of plant hormones, pollination 
and fruit setting may be subjected to 
changes due to Zn (Alloway et al., 2004). 
Our results conform with Khan et al. (2008) 
and Alam et al. (2000) who stated that the 
application of Zn enhances the 1000-grain 
weight of wheat significantly. The 
application of HA also increased 1000- 
grains weight and maximum 1000-grains 
weight was observed at 15 Kg HA ha’. The 
probable reason may be that HA increased 
the fertility of the soil which contain some 
elements which have a pivotal role in 
enhancing the fertility of the soil to a greater 
extent, thus enhancing the growth of plant 
and yield components of the crop (Khan et 
al., 2006). The increase in yield and yield 
components of wheat due to Zn application 
might be the fact that Zn played a vital role 
in the biosynthesis of the IAA and initiation 
of primodia for reproductive parts and a 
result of the favorable effect of Zn on the 
metabolic reactions within the plants Singh 
et al. (2005). 


Biological yield grain yield and harvest 
index were significantly affected by N, Zn, 
HA control vs. rest, N x Zn and N x Zn x 
HA. The treated plots showed maximum 
biological yield as compared to control. 
Biological yield, grain yield and harvest 
index were enhanced with ~ each 
augmentation of N and maximum biological 
yield was noted at 160 kg N ha’. These 
results are in contract with Khan et al. 
(2008) who stated that N_ fertilization 
improved biomass and yield of the wheat. 
Zeidan and Amany (2006) gave the same 
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results, who stated that N application 
enhances vegetative, reproductive growth 
and grain yield. Similar results are stated by 
Jan et al. (2000). Zn fertilization improved 
biological yield linearly with each increment 
of N and maximum biological yield was 
noted at 12 kg Zn ha’. The possible reason 
may be that Zn plays a vital part in plant 
growth and metabolism processes and is 
needed for enzyme activation in plants 
(Alloway et al., 2004). Our results are in line 
with Keram et al. (2014) who stated that Zn 
at the rate of 20 kg ha! combination with N 
improved equally grain and biological yield 
of wheat crop. The combined fertilization of 
N and Zn increased biological yield 
significantly. Similar results are reported by 
Asif et al. (2011) and Amanullah et al. 
(2009). Biological yield increased with HA 
increment and maximum biological yield 
was recorded at10 kg HA ha’. The possible 
reason may be that HA provides good soil 
form and enhances the accessibility of 
immobile nutrient uptake. These results 
agree with Ali et al. (2019) who reported 
that application of HA _ material has 
significantly influenced crop production, 
enhanced the fertility of soil, and minimized 
the requirement of NPK _ fertilizer 
application. 


5. CONCLUSION 


Based on the results, it is concluded that N 
fertilization at the rate of 160 kg ha’ 
produced higher grain yield and yield 
components of wheat as compared to other 
levels. Zn fertilization — significantly 
improved yield and yield contributing 
components. However, application of Zn 
above 12 kg ha! up to 18 kg ha’! was found 
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Statistically the same. HA _ significantly 
increased yield and yield components at 10 
kg HA ha’. Application of N at 160 kg ha", 
Zn at 12 kg ha! and HA at 10 kg ha" is 
recommended for obtaining higher grain 
yield and yield components of wheat. 
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ABSTRACT 

Membrane biofouling is the coverage of membrane surfaces due to 
undesirable development of biofilms causing a decrease and subsequent 
loss of productivity in water treatment settings. Continuous use of 
synthetic chemicals against biofouling is inept as it leads to the emergence 
of multi-antibiotic resistance. Application of natural products such as 
plants can be apt in curbing biofouling while checking the resistance 
challenge. This study aimedto evaluate the potential of Rosmarinus 
officinalis in the control of membrane biofouling. Bacteria from biofouling 
environments were subjected to a biofilm confirmation test and identified 
at cultural, morphological, biochemical and molecular levels. Leaves of 
R.officinalis were extracted in solvents of varying polarity and activities. 
These extracts were evaluated against bacterial biofilm formation via 
minimum biofilm inhibitory concentration (MBIC), minimum biofilm 
eradication concentration (MBEC) and mesocosm bioassays. Biofilm 
formation was confirmed in 68% of the isolates identified as 
Pseudomonas aeruginosa, Klebsiella pneumoniae and Staphylococcus 
aureus. The methanol and ethyl acetate extracts of R.officinalis indicated 
the least MICs (0.313mg/L and 1.25mg/L) against Pseudomonas 
aeruginosa and Staphylococcus aureus, respectively. Both extracts 
recorded the highest MBIC (50.00%) against Pseudomonas aeruginosa. 
The peak MBEC (57.88%) was obtained from the methanol extract 
against Staphylococcus aureus and this same extract inhibited 56.23% 
density of bacterial biofilms on glass slides. The methanol and ethyl 
acetate crude extracts of R. officinalis appreciably reduced bacterial 
biofilms; hence, this plant can be exploited as a natural antifouling agent, 
with reduced toxicity and low risk of resistance. 


KEYWORDS: Bacterial biofilms, natural product, quorum sensing, 
rosemary, water treatment. 


This is an open-access article published by the Journal of Soil, Plant and Environment (JSPAE), which 
permits use, distribution and reproduction in any medium, provided the original work is properly cited. 


1. Introduction 


Membrane Bioreactor 


www.jspae.com 


(MBR) is an 
efficient, state-of-the-art, high-quality water 
treatment technology _ that 
bioreactors modified with membrane 
filtration units for biomass retention (Meng 
et al., 2017; Waheed et al., 2017). Recent 


consists of 


technological innovations and _ significant 
footprint reduction further made MBR an 
established water treatment system (Oh and 
Lee, 2018). Despite its advantages, MBR is 
characterized with challenges such as: pre- 
treatment huddles, lack of long-term 
performance and, to a larger extent, 
membrane biofouling. Membrane biofouling 
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is the coverage of membrane surfaces due to 
undesirable development of biofilms (sessile 
slimy multicellular microbial communities) 
(Whiteley et al., 2017). This phenomenon 
hinders the effective use of MBR in water 
treatment settings (Biofilms, 2019). 


Strategies employed to control membrane 
biofouling mostly spin around physical 
cleansing of biofilms, modification of the 
membranes and incorporation of 
antimicrobial substances such as peptides 
and nitrofurazones (Hook et al., 2012). In 
addition to being pricy, the use of these 
chemicals is associated with resistance, 
environmental pollution and non-specificity 
(Lade et al., 2014; Alghamdi and Quijada, 
2019). Addressing biofoulingusing natural 
products (such as plants) as an alternative 
can be apt since the life of the bacteria must 
not be the primary target but their ability to 
form and express biofilms (Paluch et al., 
2020). Thus, it would be imperative to 
develop a strategy for mitigating membrane 
biofouling using natural products; which can 
be safe, efficient, readily available, cost- 
effective and eco-friendly. In this study, the 
leaves (commonly used part) (Meziane- 
Assane et al, 2013) of rosemary 
(Rosmarinus officinalis) were exploited 
based on its perceived ethno medicinal 
advantages and wide applications as herb 
and spice (Kalamartzis et al., 2020). The 
objectives of the research were to: (1) 
Isolate and purify bacteria from visibly 
biofouled environments, (2) screen and 
confirm the ability of the isolates to form 
biofilms and identify the biofilm-forming 
bacteria (3), extract the leaves of 
R.officinalisusing organic solvents of 
different polarity and (4) investigate the 
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antifouling activities of the extracts against 
the biofilm-forming bacteria. 


2. MATERIALS AND METHODS 
2.1. Description of Sampling Locations 


The sample sites (biofouled environments) 
were selected based on their likelihood to 
harbor biofilm-forming bacteria. These 
environments include parts of a biofouling 
model, locally fabricated at the Department 
of Environmental Science, Kaduna 
Polytechnic, Kaduna, Nigeria (10° 29' 
20.79"N, 07° 25' 21.35"E); which consist of 
membrane filter (BMm), glass (BMg) and 
plastic (BMp) substrates. Other sampling 
points include surfaces of solid objects 
(metals, plastics and wooden materials) from 
River Kaduna (RK), Kaduna, Nigeria (10° 
29' 47.13"N, 07° 25' 19.95"E) as well as 
walls and floors of water reservoirs and 
chambers (IB) in Kaduna State Water 
Treatment Plant, Malali, Kaduna, Nigeria 
(10° 33' 26.80"N, 07°29' 01.11"E). 


2.2 Sample Collection and Transportation 


A total of 117 slimy scrapings were 
collected using a simple random sampling 
technique using a sterile scoop over the 
duration of 6 months (January-June, 2020). 
Each sample was immediately transferred to 
a tightly capped Bijou bottle containing 
10mL of peptone water (Digel et al., 2018). 
The samples were transported in a cold 
condition to the microbiology laboratory, 
Department of Environmental Science, 
Kaduna Polytechnic, Kaduna, Nigeria for 
microbiological analyses. 


2.3 Isolation of Bacteria from _ the 
Biofouled Environments 
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The bacterial isolation was conducted 
using standard plate technique as described 
by Wilson et al. (2018). The biofouled 
scrapings in peptone water were vortexed 
(Digel et al., 2018) and serially diluted to 
10° using sterile normal saline as the 
diluent. A volume (100uL) from each 
dilution was inoculated (using spread plate 
technique) onto correspondingly labeled 
Nutrient Agar (gL! of peptone: 5.0, NaCl: 
5.0, beef extract: 1.5, yeast extract: 1.5 and 
agar 15; pH: 7.2) as well as MacConkey 
agar (gL of peptone: 17, proteose: 39, 
lactose monohydrate: 10g, bile salt: 1.5, 
NaCl: 5, neutral red: 0.03, crystal violate: 
0.001 and agar: 13.5; pH: 7.1) plates (Aryal, 
2019; Julistiono et al., 2018). The inoculated 
Petri dishes were incubated at 28°C for 24 
hours, following when randomly selected 
isolates were sub-cultured and purified using 
the same media, under similar culture 
conditions. 


2.4 Confirmation of Bacterial Biofilm 
Formation 


Ability of the bacterial isolates to form 
biofilms was verified using a tube method 
(TM) as described by Kirmusaoglu (2019). 
Freshly-grown bacterial culture was 
inoculated into replicate test tubes, each of 
which contained 5mL of prepared Tryptic 
Soy Broth (TSB) (gL! of tryptone: 17.0, 
soytone: 3.0, glucose: 2.5, NaCl: 5.0 and 
dipotassium phosphate: 2.5; pH 7.3 + 0.2) 
and incubated for 24 hours at 37°C. 
Following this incubation, the tubes were 
carefully emptied and the planktonic cells 
were discharged by rinsing twice with 
phosphate-buffered saline (PBS, pH 7.2). 
The sessile isolates of biofilms formed on 
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the test tubes were stained with safranin for 
1 hour. The safranin-stained tubes were 
rinsed twice with PBS to discharge the 
excess stain. After air drying, appearance of 
a visible film lining the walls and bottom of 
the tubes indicated biofilm production. The 
same volume (5mL) of a sterile uninoculated 
TSB was used as control under similar 
culture conditions. 


2.5 Identification of the Biofilm-Forming 
Bacteria 


The biofilm producing bacterial isolates 
were characterized based on _ their 
morphological properties (Gram’s _ stain 
reaction) and subjected to a series of 
biochemical (oxidase, catalase, indole, 
methyl red, Voges-Proskauer, triple sugar 
iron, citrate, urease, motility as well as H2S 
production) tests (Farinde et al., 2014; Gohet 
al., 2014; UK Standards for Microbiology 
Investigations, 2014; Cappuccino and 
Sherman, 2013; MacFaddin, 2000). 
Identities of the various bacterial species 
were confirmed using 16S rRNA gene 
sequencing (Julistiono et al., 2018) via 
Genomic DNA _ Extraction, Polymerase 
Chain Reaction, Gel Electrophoresis and 
Visualization of the PCR Products, DNA 
Band Cutting, Gel Extraction, Purification 
of the PCR Fragments and Sequencing 
(Wang et al., 2011). 


2.6 Collection and Identification of the 
Plant Sample 


Apparently healthy whole plant of 
R.officinalis was identified using standard 
keys and descriptions (Myers, 2019) and 
harvested from its cultivated site in Malali 
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Plant Gardens, Kaduna, Nigeria 
(10°32'8.7"N, 09°27'37.2" E) in May, 2020. 
A taxonomist authenticated the plant's 
identity at the herbarium section of the 
Department of Plant Biology, Bayero 
University, Kano, Nigeria. | Voucher 
specimens of the authenticated plant was 
pressed in-between clean sheets of paper, 
dried (Bulugahapitiya, 2013) and deposited 
at the herbarium for reference purpose. 


2.7 Preparation of Plant Material 


Leaves of the _ freshly collected 
R.officinalis were detached, pre-washed with 
clean water to remove extraneous materials 
(Lohaet al., 2019), rinsed with distilled 
water and distributed evenly to air-dry at 
room temperature (Gahlot et al., 2018). The 
dried leaves were excised and pulverized to 
fine powder using laboratory mortar and 
pestle. The ground powder was _ sieved 
through a 0.5mm mesh gauze to standardize 
the particle size (Teresa-May, 2018) and 
stored at room temperature in an air-tight 
dry container until needed for analyses 
(Ibrahim et al., 2017). 


2.8 Extraction of the Plant Material 


The leaf powder was extracted using cold 
maceration method (De Oliveira et al., 2019) 
in n-hexane, acetone, ethyl acetate, methanol 
and distilled water, based on the phases of 
non, less, medium, high and comparably 
high polar solvents, respectively 
(Bulugahapitiya, 2013). According to a 
procedure described by Sagbo et al. (2020), 
100g of the leaf powder was macerated in 
1000mL of the solvent and the set-up was 
allowed to stand for 72 hours at room 
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temperature with intermittent agitation (De 
Oliveira et al., 2019; Loha et al., 2019). The 
damp plant material was passed through a 
cheese cloth, allowed to settle and re-filtered 
via a Whatman grade | filter paper (1 1um) 
(Loha et al., 2019). The residue (marc) was 
re-extracted in similar solvent to recover as 
much occluded solution as possible. The 
extracts were concentrated using rotary 
evaporator at 40°C in vacuo and air dried in 
a fume hood (Teresa-May, 2018). The dried 
fractions of the crude extracts were stored in 
air tight glass containers at 4°C under 
refrigeration until required for further 
analyses (Sagbo et al., 2020). 


2.9 Preparation of the Plant Extracts 


Exactly 0.1g of the dried crude extract was 
dissolved in 10mL of 1% Dimethyl 
Sulfoxide (DMSO) to achieve 10mg/mL 
(Famuyide et al., 2019). This (stock) was 
diluted serially to prepare subsequent test 
concentrations (5.000, 2.500, 1.250, 0.625, 
0.313, 0.156 and 0.078mg/mL). Sterility of 
the extracts was verified by inoculation (via 
streaking) on freshly prepared NA, which 
was incubated at 37°C for 24 hours. 


2.10 Minimum Inhibitory Concentration 
(MIC) Assay 


Five milliliters (SmL) of standardized (0.5 
ODs9snm) bacterial inoculums (Lade et al., 
2014) from an overnight culture was 
inoculated into 95mL of freshly prepared LB 
broth (Lade et al., 2014) and incubated at 
37°C for 48 hours (Julistiono et al., 2018). 
Wells of a flat bottom polystyrene 
microtiter plate were conditioned by 
introducing 200uL of (plain) LB broth to 
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each and allowed to stand for 1 hour at room 
temperature (Lade et al., 2014). The wells 
were emptied and 20uL of the activated 
culture, followed by 180uL of the plant 
extract (0.078—10.000mg/mL)were added 
to each (Taufik et al., 2018). The plate 
was covered and incubated under static 
condition at 28°C for 48 hours (Biswa and 
Doble, 2013). Plain (uninoculated) LB 
broth was used as negative control while 
ciprofloxacin (5ug/mL) was applied as 
positive control. Bacterial MIC was 
determined with the use of a micro plate 
reader (AMR-100, China) at ODs59s as the 
minimum concentration where 
absorbance of the treatment was less or 
equal to that of the negative control (Da 
Rosa et al., 2016). 


2.11 Minimum Biofilm Inhibitory 
Concentration (MBIC) Assay 


Inhibition of the bacterial biofilm 
formation was determined using a crystal 
violet method in microtiter plate (Julistiono 
et al., 2018). Equal volumes (100uL each) of 
the activated culture and the extract 
(0.078-1.250mg/mL) were introduced to 
each well following conditioning. The plate 
was incubated in static condition at 28°C for 
48 hours. Supernatants from the wells were 
carefully aspirated out without disrupting 
the biofilms on the base and the wells 
were washed thrice with PBS to remove 
any unattached bacterial cells. The plate was 
incubated at 37°C for 15minutes. Formed 
biofilms were fixed with 200uL of 99% 
methanol for 20 minutes and stained with 
100uL of 0.2% (w/v) crystal violet solution 
for 15 minutes at room temperature. Excess 
stain was removed from the wells by 
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rinsing four times with PBS, after which 
100uL of 95% ethanol was introduced to 
extract the crystal violet in solution from 
the biofilms (Hossain et al., 2017). 
Absorbance of the dissolved crystal 
violet (which corresponds to a measure 
of bacterial cells that formed the 
biofilms) was determined at 595nm 
(Taufik et al., 2018) using the micro 
plate Reader. Inoculated LB broth 
(without the extract) was used as control. 
The MBIC was determined as_ the 
minimum concentration where 
absorbance of the treatments was less 
than or equal to that of the control (Da 
Rosa et al., 2016). Percentage inhibition 
of biofilm was calculated using the 
equation below (Julistiono et al., 2018): 


% biofilm inhibition (595nm) 
= OD of control — OD of test 


OD of control oy 


2.12 Minimum Biofilm Eradicating 
Concentration (MBEC) Assay 


Following conditioning of the micro 
plate wells, 100uL of the activated culture 
was inoculated into each well. The set-up 
was incubated under static condition at 28°C 
for 48 hours (Taufik et al., 2018). The 
liquid culture was aspirated out and the 
wells were washed with PBS. Wells with 
successfully induced biofilms were filled 
with 200uL of the extracts at 0.078-— 
1.250mg/land the plates were incubated 
for 6 hours at 37°C. The biofilms were 
fixed with 200uL of 99% methanol for 
20 minutes and stained with 100uL of 0.2% 
(w/v) crystal violet solution for 15 minutes 
at room temperature. Excess stains were 
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removed by rinsing four times with PBS, 
which preceded the addition of 100uL of 
95% ethanol (to each well) to extract the 
crystal violet in solution from the biofilm 
(Hossain et al, 2017). Absorbance 
(595nm) of the dissolved crystal violet 
was determined with the use of micro 
plate reader as a measure of cells that 
formed the biofilms (Taufik et al., 2018). 
Inoculated LB broth (without the extract) 
was used as control. The MBEC was 
calculated as the minimum concentration 
where absorbance of treatment was at 
least 50% less than that of the control, 
indicating up to 50% eradication of the 
formed biofilms (Da Rosa et al., 2016). 


2.13 Mesocosm Experiment 


To further confirm their antifouling 
activities, the crude plant extracts were 
subjected to mesocosm experiment in 
accordance with a procedure described by 
Dobretsov et al. (2011). Specifically, 2 sets 
of 2L capacity transparent plastic containers 
were filled with IL of the compound 
solution. This solution was prepared by 
dissolving the extract in standard unfiltered 
seawater, sampled from the Eleko beach, 
Lagos, Nigeria (06° 26' 17.23"N, 03° 51' 
06.51"E) using a_ standard procedure 
(Cheesbrough, 2006). Final concentrations 
of 0.313mg/mL and 0.156mg/mL extract 
in seawater were used. Replicates of sterile 
microscope glass slides (25 x 75mm) were 
immersed horizontally into each container. 
Microscope slides were also dipped into a 
similar container filled with an equal volume 
of the unfiltered seawater (control). 
Containers with the slides were kept under 
illumination at room temperature 
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(25°C+2°C) for 5 days (Wilson et al., 2018). 
The slides were brought out and fouling was 
fixed with 1% formaldehyde in seawater. 
These slides were then stained with DNA- 
binding fluorochrome  4,6-diamidino-2- 
phenylindole solution (0.5ugmL"') and air 
dried. Bacteria in 10 randomly selected 
fields of view were enumerated under an 
epifluorescence microscope (Wild M20, 
Switzerland) using direct count at a total 
magnification of x2000. 


Percentage biofilm inhibition 
control count — test count 


= —— /> 100 
control count 


2.14 Statistical Analyses 


Data were presented as mean + standard 
deviation (SD) of replicate assays. The 
mean and standard deviation of bioassays 
were computed using Microsoft Excel 
(version 2016). Values of inhibitory 
activities were appraised by one-way 
analysis of variance (ANOVA) with the 
use of GraphPad Instat (version 3.10); in 
comparison with controls. All P-values 
<0.05 were regarded as _ statistically 
significant, which were illustrated with 
different superscript alphabets, while those 
>0.05 were considered insignificant and 
denoted by similar superscripts. 


3. RESULTS AND DISCUSSION 


In this study, biofilm production among 
bacteria isolated from the various biofouled 
environments was rated ‘strong’ where there 
was a visibly high biofilm adherence to both 
the wall and bottom of the test tubes. It was 
regarded as ‘moderate’ in case of less 
adherence to the tubes and ‘weak’ where 
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only a trace of biofilms manifests in the 
walls or bottom of the test tube(s). The 
result indicates that the trait of biofilm 
formation was confirmed in the majority 
(63.16%) of the screened isolates. This 
implies that most of the sampled 
environments harboured the typical biofilm- 
forming bacteria, which can easily be 
isolated using protocols adopted in this 
research. Among the biofilm formers, 
isolates from the membrane filter, plastic 
and glass substrates (Mi, M2 and Ms3, 
respectively) of the biofouling model 
excelled with up to 25% ‘strong’ formation. 
Since all the three isolates originated from 
the (controlled) biofouling model, this 
signifies that biofilm formation could 
environmental 
conditions such as availability of nutrients, 
moisture and temperature. The result of 
Zuberi and Nadeem (2017) corroborates the 
finding of this study as they similarly 
reported as high as 63.64% biofilm 
formation in bacteria isolated from contact 
lenses and their accessories in Karachi, 
Pakistan. 


depend strongly —=— on 


4.1 Identity of the Biofilm-Producing 
Bacteria 


The isolates with confirmed biofilm 
formation ability were identified as 
Staphylococcus aureus RBSB2_Cl, 
Pseudomonas aeruginosa NT 10038 and 
Klebsiella. pneumoniae (C2244. This 
confirmed the inherent biofilm formation 
property of these bacterial species. Awoke et 
al. (2019) support this finding as they 
similarly identified the capability of biofilm 
production in Staphylococcus aureus, 
Pseudomonas aeruginosa and Klebsiella 
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sp. in Southwest Ethiopia. Although, 
Staphylococcusaureus itself has been 
reported (Paluch et al., 2020) to control 
biofilm production in some other bacteria. 


4.2Minimum Inhibitory Concentration 
of the Plant Extracts 


The MICs of the extracts were 
determined as 0.31 to 1.25mg/mL against 
Klebsiella 
pneumonia and Staphylococcus aureus. 
Majorly, lower MICs were recorded from 
the methanol and ethyl acetate extracts. 
This might owe to the fact that many 


Pseudomonas aeruginosa, 


phytochemicals are polar in nature; 
hence, significant biological activities 
concentrate mostly in the polar regions 
(Bogavaca et al., 2017). The least MIC 
(0.31mg/mL) was obtained from the 
methanol extract against Pseudomonas 
aeruginosa, followed by 0.63mg/mL 
from the same extract against Klebsiella 
pneumonia. This implies that the biofilm 
producing bacteria were most susceptible 
to the methanol extract. The study of 
Van-Vuuren (2008) corroborates this as 
MICs <1.25mg/mL of some South African 
plants were reported and regarded as strong 
values. In the result of Bogavaca et al. 
(2017), much higher (5Omg/mL) MIC value 
of R. officinalis (but essential oil)was 
identified against Pseudomonas 
aeruginosa. Similar to our finding, Jarrar 
et al. (2010) recorded activities of the 
ethanol extract of R. officinalis, collected 
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Table 1: Confirmed Biofilm Production of Bacterial Isolates from Slimy Surfaces of the Membrane Filter, Plastic and Glass Substrates of the 
Biofouling Model, Water Reservoirs from Kaduna State Water Treatment Plant and Solid Objects from River Kaduna 


Bacterial Isolates kb Gb dl Is Ie Ib Mi M2 Ms Ma Ms Ri Ro R3 Ra Rs Ro Ro 
Extent of Biofilm Production - ++ ++ + - + ++ +++ 44+ 444+ - ++ 4+ - - + + = = 
Percentage Biofilm Formation: N (%) 19(63.16) 
Strong: n (%) 12(25.00) 
Moderate: n (%) 12(33.33) 
Weak: n (%) 12(41.67) 
Key: +, ++ and +++ = Weak, moderate and strong biofilm production respectively; - = No visible biofilm production. 


M:-s = Isolates from plastic, glass and membrane filter substrates of the biofouling model, R1-7 = Isolates from biofouled substrates from River 
Kaduna and [1-7 = Isolates from surfaces of water reservoirs from Kaduna state water treatment plant. 
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Table 2: Cultural, Morphological and Biochemical Identities of the Biofilm-Forming Bacteria 


Colonial Characteristics Gram’s Ox Cat Ind MR VP Gl Lc Su H2S Cit Ur Mot Inference 
Stain 

on NA on MacConkey Reaction 

Blue-green Yellow-green Gram - ++ +++ - - - + + - - + - + Pseudomonas 
bacilli aeruginosa ** 

Mucoid milky Pink mucoid Gram - - + - - a i - + + -  Klebsiellasp. 
bacilli 

Large smooth NG Gram + - + - + + + + + - + + - Staphylococcus 

circular cocci aureus 


Key: - = negative, + = positive, Ox = oxidase, Cat = catalase, Ind = indole, MR = methyl red, VP = VogesProskauer, Gl = glucose, Le = 
lactose, Su = sucrose, Cit = citrate (Simon’s), Ur = urease, Mot = motility, ** = grown at an elevated temperature (42°C), NG = No growth. 


Table 3: Minimum Inhibitory Concentration (mg/mL) of Rosmarinus officinalis Leaf Extracts against Pseudomonas aeruginosa, Klebsiella 
pneumoniae and Staphylococcus aureus 


Pseudomonas aeruginosa Klebsiella pneumoniae Staphylococcus aureus 
Extracts Aq M E A H Aq M E A H Aq M E A H 
R.officinalis 2.500 0.313* 1.250 - - - 0).625* - - 2.500 - - 1.250 - - 
Ciprofloxacin < 0.078 < 0.078 < 0.078 


Plain LB = = = 


Key: Aq = Aqueous, M = Methanol, E = Ethyl acetate, A = Acetone and H = N-hexane. Results <1.250mg/mlL (especially those indicated with*) 
were considered strong MIC values. - = MIC values >5.000mg/mL.Ciprofloxacin and Plain LB = positive and negative controls respectively. 
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Table 4: Inhibitory Potential of the Methanol and Ethyl Acetate Leaf Extracts of Rosmarinus 
officinalis Against the Biofilms of Pseudomonas aeruginosa,Klebsiella pneumonia and 
Staphylococcus aureus 


Inhibition of Biofilm Development (%) 


Bacterial species 


Extracts Pseudomonas Klebsiella Staphylococcus aureus 
(0.078—1.250mg/mL) aeruginosa pneumonia 

Methanol 50.00+0.66(b)® 43.10+3.84(c)! 37.69+5.48(c)! 

Ethyl Acetate NA 42.414+2.24(a)! 48.8+1.00(b)! 
Ciprofloxacin 51.90+12.61" 51.3843.51' 53.8044.82* 


(0.078mg/mL) 


Values are mean (+SD) percentage biofilm inhibition. 

Values with different superscripts across the same column are significantly different (P<0.05). 
Key: NA = not active against biofilms at all tested concentrations. 

a, b, c, d and e = 0.078, 0.156, 0.313, 0.625 and 1.250mg/mL respectively. 


Table 5: Eradication Potential of Rosmarinus officinalisLeaf Extractsagainstthe Biofilms of 
Pseudomonas aeruginosa,Klebsiella pneumonia and Staphylococcus aureus 


Eradication of formed Biofilms (%) 
Bacterial species 


Extracts Pseudomonas Klebsiella Staphylococcus 
(0.078—1.250mg/mL) aeruginosa pneumoniae aureus 
Methanol 48.67+6.14(c)™ S11 7t2.24(6)" **57.88+4.76(c)™ 
Ethyl Acetate FS) SOLO. 7 9(C)" 41.92+4.09(d)" 48.22+3.84(c)" 
Ciprofloxacin **64_21+43.02? **61 5543.38? **56.3941.58? 
(0.078mg/mL) 


Values are mean (+SD) of percentage biofilm eradication. 

Values with different superscripts across the same column are significantly different (P<0.05). 
Key: - = no biofilm eradication recorded, ** = active eradication of biofilms (=50%) 

a, b, c, d and e = 0.078, 0.156, 0.313, 0.625 and 1.250mg/mL respectively. 
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Table 6: Density of Biofilm-Forming Bacteria on Glass Slides Exposed to the Extracts of 


Rosmarinus officinalisin Unfiltered Seawater 


Count (cells/mm’) 


Plant Extracts Control 0.156mg/mL 0.313mg/mL 
Methanol 127.93+4.654 102.47+2.96' (19.90%) 56.00+2.76* (56.23%) 
Ethyl Acetate 125.67+4.194 103.17+3.59' (17.90%) 73.00+4.16° (41.91%) 


Values are mean (+SD) bacterial density with percentage inhibition in parenthesis. 
Mean values with different superscripts across the same row are significantly different (P<0.05). 


from the Northern Palestine against 
Staphylococcus aureus, at an MIC of 
0.39mg/mL. It was established (Abkhoo et 
al, 2010) that ethanol extracts of R. 
officinalis from Tehran, Iran, significantly 
inhibited the growth of Pseudomonas 
aeruginosa at the MIC of 0.1mg/mL. 


Generally, variation in MIC values was 
observed across the three (3) bacterial 
isolates in the present study. This can be due 
to the presence of different intrinsic levels of 
tolerance to the tested plant compounds, as 
similarly observed by Ahmad and Aqgil 
(2007). 


4.3 Biofilm Inhibitory Potential of the 
Plant Extracts 


Almost all the extracts imbhibited the 
formation and development of preformed 
biofilms at appreciable limits. The extracts 
were able to hinder the development of 
biofilms formed by Pseudomonas 
aeruginosa, Klebsiella pneumonia and 
Staphylococcus aureus at the range of 24.89 
to 50.00%. The methanol extract displayed 
the highest (50.00%) potential of biofilm 
inhibition against Pseudomonas 
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aeruginosa at 0.16mg/mL. This unveils the 
puissance of R. officinalis as a potential 
antifouling agent. The probable mechanism 
of this biofilm inhibition may be reduction 
in the production of extra polymeric 
substances (EPS), which is an important 
component of bacterial biofilms, crucial to 
the maintenance of the spatial structure of 
the consortium, as established by Paluchet 
al. (2020). Exactly 50.00% of biofilm 
formed by Pseudomonas aeruginosa was 
similarly inhibited but at 7.80mg/mL of R. 
officinalis in the study of Yazdeliet al. 
(2021). Ziemichéd and Skotarczak (2017) 
further established the capability of plant 
products to inhibit both formed and 
preformed biofilms of Staphylococcus 
aureus. Studies conducted by Endo et al. 
(2018) were in agreement with ours as they 
equally revealed as high as 50.00% activity 
against the preformed biofilm = of 
Staphylococcus aureus and Pseudomonas 
aeruginosa at 30.00 to 250.00ug/mL of the 
leaf extracts of R. officinalis. Likewise, 
Ceylan et al. (2014) recorded 39.49% and 
51.30% inhibition capacity of the preformed 
biofilms of Pseudomonas aeruginosa and 
Staphylococcus respectively 
using R. officinalis essential oil at 10.00 to 


aureus, 
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0.08ug/mL. Up to 57.00% of biofilm formed 
by Staphylococcus  Epidermidis was 
inhibited by the essential oil of R. officinalis 
obtained from Tunisia, as reported in the 
work of Jardak et al. (2017). 


4.4 Biofilm Eradication Potential of 
the Plant Extracts 


Results of the eradication assay 
revealed that biofilms formed by 
Klebsiella 
pneumoniae and Staphylococcus aureus 
were appreciably reduced by the extracts, 
indicating an activity range of 33.95 to 
57.88%. The highest (57.88%) activity was 
obtained from the methanol extract of R. 
Officinalis at O0.31mg/mL against the 


Pseudomonas aeruginosa, 


biofilm formed by Staphylococcus aureus. 
Jardak et al. (2017) also established that the 
biofilm formed by — Staphylococcus 
epidermidis was eradicated by up to 67.53% 
when exposed to R. officinalis essential oil 
at the concentration of 50.00mg/mL. In this 
study, R. officinalis might have reduced 
bacterial biofilms by producing and 
releasing compounds capable of affecting 
bacterial molecular signals and inhibiting 
behaviors under the control of quorum 
sensing (QS) (Yazdeli et al., 2021). 


4.5 Effect of the Extracts on the 
Density of Biofilm-Forming Bacteria 


From the mesocosm experiment, both 
methanol and ethyl acetate extracts of R. 
officinalis inhibited the formation of 
microbial communities on the glass slides at 
both tested concentrations. The methanol 
extract significantly (P<0.05) decreased 
bacterial densities (56.23% and 19.90%) 
at 0.313mg/mL and  0.156mg/mL 


www.jspae.com 


Salisu et al. 


respectively; in relation to the control. 
Activities recorded from this extract 
(especially at 0.33mg/mL) might indicate 
probable inhibition of QS among the 
exposed bacteria. This might have led to low 
attachment, hindering the subsequent 
biofilm formation and development. This 
idea was supported by Dobretsovet al. 
(2007), who established that QS inhibitors 
affect microbial composition and densities. 
According to Kjelleberg et al. (2001), anti- 
quorum sensing (AQS) agents may alter 
bacterial composition in biofilm formation, 
leading to a shift in the microbial 
communities from being dominated by 
Gram-negative to Gram-positive bacteria. 
However, this does not exclude the 
possibility that the biofilm inhibition 
(recorded in this study) might have also 
arisen from the obstruction (by the extracts) 
of other regulatory cascades that may govern 
the process of biofilm formation, as 
highlighted in the study of Dobretsovet al. 
(2011). 


5. CONCLUSION 


The biofilm formation trait was confirmed 
in the majority of the bacteria isolated from 
the biofouled environments. These bacteria 
were identified as Staphylococcus aureus 
RBSB2_Cl, Pseudomonas aeruginosa 
NT10038 and Klebsiella pneumoniaeC2244. 
Evaluation of the antibiofilm activities of R. 
officinalis revealed that this plant (especially 
its methanol extract) possessed some 
features of inhibitory activities against 
bacterial biofilms. This indicates that the 
methanol, followed by ethyl acetate extracts 
of this plant, might contain many 
phytoconstituents at varying concentrations, 
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which can be responsible for their 
antifouling effects. Therefore, this study 
has identified the potential of R. 
officinalis as, a source of active 
antifouling compounds that can be cheap, 
eco-friendly and readily available. 
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Research Farm, Department of Agronomy, University of Agriculture 
Faisalabad, Pakistan, in 2015. The experiment consists of three sorghum 
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13.66), stem diameter (1.09 cm), forage yield (55.52 ton ha"), dry matter 
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Pl and P2. This study suggested that the cultivar sorghum 2011 with a 
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1. INTRODUCTION important agriculture sector and has a 
central role in our country's rural economy. 
Its accounts for 54% of agricultural GDP 
and 11.9% of total GDP (Mehmood et al., 
2020, Jahnke and Jahnke, 1982). Sorghum 
(Sorghum bicolor L.) is a member of the 
Poaceae family is locally called Jawar and is 
commonly used as fresh, silage or hay form. 


Agriculture is the most important sector 
of Pakistan and greatly influences economic 
growth. It accounts for 21.4% of GDP and 
45% of employment is engaged with 
agriculture (Gecho, 2017, Iqbal et al., 2015). 
Development in the agriculture sector 
stimulates growth in the agro-industry, 


j : ; ; Sorghum is a dual-purpose (i.e., grown for 
especially the textile sector. Livestock is an e burp ( e 


grain and fodder yield) summer season crop 
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(Hedayetullah and Zaman, 2018; Sarfraz et 
al., 2012). It requires fewer resources and 
gives a high yield being the most drought- 
tolerant crop of the universe (Ringo et al., 
2014) and can be grown in tropical and sub- 
tropical countries of the world. Sorghum has 
more potential to fulfill the future demands 
of livestock, grain food and beverages. Its 
syrup extracted from sweet sorghum is 
mainly used in ethanol production (Klasson 
et al., 2021, Rao et al., 2016; Kirouani et al., 
2021). Sorghum is the fourth cereal crop of 
Kharif season and an important forage crop 
in many regions of the world, including 
Pakistan (Ghani et al., 2015; Rana et al., 
2014). It exhibits rapid growth, is relatively 
resistant to dryness, high productivity and 
high percentage of protein. In short period 
sorghum produced a large amount of seed 
and fodder (Djanaguiraman et al., 2020). Its 
fodder comprises 11 % protein, 71 % 
carbohydrates, 2% crude fiber, mineral and 
nitrogen-free extract (Hussain et al., 2020). 
Sorghum accounts for half of the forage in 
the rainfed area. In Pakistan sorghum was 
cultivated on 198thousand hectares land 
with the production of 123 thousand tones 
and an average grain yield of 621 kg ha 
‘(Hussain et al., 2015). 


Fodder crops play a crucial role in the 
agricultural economy of developing 
countries by providing the cheapest source 
of feed for livestock (Herrero et al., 2013, 
Upton, 2004; Ramana, 2022). Livestock is a 
vital part of farming plays an important role 
in the economic development of the rural 
community of Pakistan. Livestock accounts 
for 55 % of agricultural GDP and 12 % of 
total GDP (Hussain et al., 2015). Lower 
fodder production and less accessibility to 
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feed are the main factors of decreasing 
livestock in Pakistan. Furthermore, 
providing quality animal feed in a suitable 
amount can increase livestock production. 
Fodder production fulfills 30 to 50 % of the 
total fodder consumption in Pakistan. While, 
low quality of animal feed caused low meat 
and milk production (Nouman et al., 2014). 


The available fodder contributes 1/3 
less than that needed fodder and _ its 
deficiency is further increased due to a 
decline in area under fodder crops by 2% 
after each decade (Nadeem et al., 2017, 
Herrero et al., 2013). Amongst the Kharif 
forage crops, sorghum is an important one 
that possesses a wide range of ecological 
flexibility. The extension growers largely 
sow it for feed and fodder in rainfed and 
irrigated regions of the country. Almost 
sorghum is fed to every class of livestock 
(Kumar and Upadhyay, 2008). The 
performance of dairy animals depends on 
the regular availability of quality fodder 
insufficient amount. The vital limitation on 
profitable animal production in developing 
countries is the unavailability of quality 
forage (Kumar and Upadhyay, 2008, 
Quddus, 2012). Considerable, differences 
have been reported among the sorghum 
cultivars for yield, quality traits (Rao et al., 
2013), and response to planting densities 
(Moosavi and Sciences, 2012). 


Keeping in view the importance of 
fodder crop, the present study was 
conducted to find out the most suitable 
sorghum variety under suitable planting 
pattern for higher yield and fodder in terms 
of quality and quantity in the agro- 
ecological climate of Faisalabad. 
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2. MATERIAL AND METHODS 
2.1 Experimental Site and Design 


The experiment was performed at 
Agronomic Research Area, University of 
Agriculture Faisalabad, Pakistan in 2015. 
The experimental site lies between 30.35- 
41.47°N latitude and 72.08-73.40°E 
longitude at an elevation of 184.4 m above 
sea level. The experimental site's maximum 
and minimum temperature, rainfall, sunshine 
hours, relative humidity, wind speed, and 
evapo-transpiration were recorded at local 
meteorological _ station. The weekly 
maximum and minimum _ values of 
temperature, rainfall, sunshine hours, and 
wind speed are given below in figure 1. 

The experiment was laid out in 
randomized complete block design with split 
plot arrangement using three replications 
and a net plot size of 6 m x 7.2 m. All the 
varieties were sown with a seed rate of 75 
kg ha’. Fertilizers are applied at the rates of 
58:58:0 (N: P: K). The experiment 
wasconsisted of three Sorghum cultivars 
(Jawar 2002, Sorghum-2011, and JS-2002) 
and three treatments of planting patterns 
(P1l= 60 cm x 20 cm, P2 = 50 cm x 24 cm 
and P3= 40 cm x 30 cm). All other 
agronomic practices were kept the same for 
all the treatments. 


2.2. Measurement and Analysis 
2.2.1. Agronomic and _ Yield-related 
Parameters 


Sorghum plants were counted in the one- 
meter length of three randomly selected 
rows in each plot and then averaged per 
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square meter. For plant height (cm), ten 
sorghum plants were selected and_ their 
height was measured from the base to the tip 
of the longest leaf with measuring tape and 
then averaged. The total number of leaves 
from ten plants was counted and then 
average leaves per plant were calculated. 
Furthermore, to determine leaf area per plant 
(cm’), at each harvest, leaves were removed 
from ten randomly selected plants and 
passed through the leaf area meter model LI- 
3000 and readings were noted. For the 
measurement of stem diameter (cm), the 
diameter of ten randomly selected plants 
from each plot was measured with the help 
of Vernier Caliper from the base, middle and 
top portions of the stem and then averaged. 


Moreover, for the determination of the 
weight per plant (g), ten plants were 
randomly selected from each plot at each 
harvest with the help of sickle. Each plant 
was weighed and averages of these plants 
weights were calculated to get the fresh 
weight of each plant in grams. Fresh weight 
per plant (g) was observed by selecting five 
plants randomly and taken from each plot 
then weighted to determine the mean fresh 
weight per plant. While to evaluate dry 
weight per plant (g), fresh samples were 
dried at 60°C for 48 hours in a fan-assisted 
oven until a constant weight was reached 
and weighted to obtain the mean dry weight 
per plant. For the determination of forage 
yield (t ha'!), all the crop plants in each net 
plot reserved for recording yield at final 
harvest and weighed separately with the help 
of a spring balance and converted into t ha”. 


91 


Journal of Soil, Plant and Environment 


2 
@ aie 
= 
e 15-4 
® 
a 
= 
~~ 10 
5-34 
Oo 7 - y ’ ’ Y 
a ot at at at at gt gt at st st st st 
gS FFF 
PAS?’ HP ps GE Ee aF gF GF gGS_ SF as 


Asim et al. 


6.0 
5.5 
rN 
5% ‘ / 
ry - 2 
: / 50 § 
/ a 
i \ f = 
j nN se) 
; | 45 £ 
' \ d J _ 
\ \;/ “” 
é f y = 
\ / 40 8 
1 
, f 
\/ 3.5 
' ' ' ' —_-o-, $- 3.0 
o oF oF oF ot of ot of of of of of of 
Fe 
PagP FF FF GF aS SF BESS AS AS aS 


Weeks after sowing 


—— Weeks after sowing vs Max T 
Weeks after sowing vs Min T 


( J Rain fall 
Weeks after sowing vs SSH 
Weeks after sowing vs Vind 


Figure 1. Weekly maximum and minimum values of temperature, rainfall, sunshine hours, and 


wind speed of the experimental site 

The dry matter yield (t ha!) were assessed 
by selecting randomly ten plants at harvest 
from each plot and chopped with the help of 
a forage cutter and then thoroughly mixed. 
The fresh weight of the sample was recorded 
and a sample of 500g was taken from each 
plot and dried in an oven at 70°C to a 
constant dry weight. These plants were 
selected from the plot area used for green 
forage yield and their weight was added in 
each respective plot. Dry matter percentage 
calculated was used to convert green forage 
yield to dry matter yield. Furthermore, to 
record dry matter (%), a chopped known 
weight of forage from each plot was taken 


and then dried at 80 °C in an electric oven to 
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a constant weight. The dry matter 
percentage for each plot was calculated 
following the formula below. 

Dry matter (%) = Dry weight/fresh 
weight x 100 
2.2.2. Quality Parameters 
3.2.1. Crude Protein (%) 

Initially, the samples were grinded with a 
locally made grinder until the sample was 
almost converted into a powdered form and 
no sieve was used. The powdered sample 
(1g) was added to KJeldahl digestion flask 
along with 30ml concentrated H2SO4 and 
10g digestion mixture. The powdered 
sample (1g) was added to KJeldahl digestion 


flask along with 30ml concentrated H2SO4 
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and 10g digestion mixture. After keeping for 
half an hour it was heated slowly in the 
beginning and then on full heat until a 
transparent green liquid material resulted. 
On cooling transferred to a 250 ml 
volumetric flask and volume was made up to 
the mark. Aliquot of 10 ml from this 
material was taken in the micro kjeldahl 
apparatus using 15 ml 40% NaOH for each 
sample. It involves digestion of the plant 
material dried at 70°C with concentrated 
sulphuric acid and digestion mixture, 
comprising K2SO4, CuSO, and FeSO, in the 
ratio of 10: 0.5:1. 

Nitrogen evolved as ammonia was 
collected in a receiver containing boric acid 
(4%) solution and mixed indicator of 
bromocresol green and methyl. The 
distillation was titrated against N/10 
sulphuric acid till the original color of 
methyl red was restored. Blank was run to 
eliminate the percentage of nitrogen present 
in other chemicals used to digest the sample. 
From the quantity of acid used in titration, 
the percentage of element nitrogen was 
calculated by using the formulas. The 
reading obtained was multiplied by 6.25 to 
get crude protein percentage. The crude 
protein percentage was determined by using 
the standard procedure as recommended by 


(Salo-vaéaéndnen and Koivistoinen, 1996). 
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N (%)=A-Bx100x100x0.0014/ Volume of 
digested sample used. 

Where A= quantity of acid (N/10 H2SOu) 
used. 

B= Blank reading (N/10 H2SO4 used in 
blank reading), 100= volume made after 
digestion, 100 for percentage (Which is 
equal to grams of N in Iml of N), 

0.0014= Factor (Which is equal to grams of 
N in Iml of N/ 10H2SO4 

2.2.3. Crude Fiber (%) 

Two grams of oven-dried sample were 
digested in 200 mL of 1.25 % H2SO4 in 500 
ml beaker for 30 minutes to determine crude 
fiber. Then contents were filtered by linen 
cloth and residues were washed and digested 
again with 200 ml 1.25% NaOH for 30 
minutes and after that, it was again filtered 
and washed. The residues were put in a 
weighed china dish and dried in an oven for 
24 hours at 105°C. After recording, the dry 
weight samples were placed in a muffle 
furnace at 600°C until grey or white ash was 
obtained. The weight of the ash was 
recorded. Crude fiber (%)= {(Wt. of dried 
residues — Wt. of ash)/Wt. of the dried 
sample)} x100 

The crude fiber percentage was 
determined by using the standard procedure 
as recommended by (Salo-véanénen and 


Koivistoinen, 1996). 
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2.2.4. Ash Contents (%) 

A 5g of oven-dried sample was placed in 
a clean previously weighed china dish (W1) 
to determine ash content. The samples were 
placed in a muffle furnace at (550-650°C) 
until white or grey ash was obtained. After 
that, residues were cooled in a desiccator 
and recorded the weight (W2), and the 
percentage was calculated as follows: 
Total ash % = [(W2 — W1) + (Weight of the 
sample)] x 100 

Total ash percentage was determined 
using the standard procedure proposed by 
AOAC (1990). 
2.3. Statistical Analysis 

Data collected on all parameters were 
analyzed _ statistically using MSTAT-C 
software (Crop and = Soil Sciences 
Department of Michigan University of the 
United States). The least significance 
difference (LSD) test at the 5% probability 
level was applied to compare the treatment's 


means (Shrestha, 2019). 


3. RESULTS 
3.1 Growth and Yield Traits 

Our results showed that growth and yield 
attributes were significantly affected by 
different different 


planting _ patterns, 
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cultivars and their interactions (Table 1). 
Higher planting density (40.67 m7), plant 
height (237.11 cm), leaves per plant (13.67), 
leaf area per plant (2520.20 cm’), stem 
diameter (1.11cm), weight per plant (280.94 
g) were recorded in Sorghum 201 |cultivar 
as compared to other cultivars. In the case of 
planting patterns, P3 (40 x 30) resulted in 
higher planting density (40.77 m”), plant 
height (236.33 cm), number of leaves 
(13.66), stem diameter (1.09 cm), weight per 
plant (136.3 g) and dry weight per plant 
(40.34g) as compared to P; and P2.The 
lowest growth and yield attributes were 
recorded in P;(60 x 20 cm). 

Data on forage yield showed significant 
difference among cultivars and planting 
patterns (Figure 2). The cultivar Sorghum- 
2011 produced maximum forage yield (57.6t 
ha!) followed by Jawar-2002 (52.5 t ha‘') 
and cultivar JS-2002 produced a minimum 
forage yield (45.5 t ha"). In case of planting 
patterns, forage yield were found significant 
and it ranged from 47.39-55.5 t ha’!. The 
maximum forage yield (55.52 t ha) was 
found in plots where sorghum was sown 
using narrow row spacing P3 (40 x 30). The 
minimum forage yield (47.39 t ha!) was 


recorded in P; followed by P2. 
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Table 1. Response of Sorghum planting density (m7), plant height (cm), leaves plant’, leaf area 


plant'stem diameter (cm), weight plant (g), fresh weight per plant (g) and dry weight per plant 


(g) to different varieties and planting patterns. 


Fresh Dry 
Planting Nad Pane Leaves LA Stem Weight weight weene 
Varieties | density PH (cm) 1 1 diameter per 
Pattern 2 plant plant plant (g) per 
(m™) (cm) plant 
plant (g) 
(g) 
60x20 ) Jawar | 3933h 229.67cd 12.00b 2224.80 1.04ab 241.45d 132.74c 38.00bc 
cm 2002 
ae ae 39.00b 234.00c 12.33b 2406.50c 1.06ab 251.23c_ ~—-135.03be 38.89be 
60 x 20 
cm | 38-2002 | 37.33¢ —223.67d_~—11.00e 2206.00c 1.08ab_—-235.36d—:138.22b 35.51c 
50x24) Jawar | 39435 232.67¢ 13.00ab 2406.50c 1.09ab -261.23b_135.03be 38.89bc 
cm 2002 
aang aaa 41.00a 236.33b 13.67ab 2513.40b 1.10a 265.14b 138.68  40.40ab 
50 x 24 
cn | SS-2002 | 38.33b  229.33cd 11.66¢ 2240.40e —1.13a_252.23e 141.592 42.53a 
40x30 | Jawar | yi oa 235.00b14.00a_2470.90c_ 1.02ab -283.45ab-—«138.22b 40.08ab 
cm 2002 
ied a 42.00a 241.00a 15.00a 2640.70a 1.02ab 294.14a 141.592 42.53a 
40 x 30 
cm | §S-2002 | 39.33b —233.00c_—:12.00b 2340.30d_-1.04a—-265.24b 129.08c _38.40bc 
ANOVA (CF values) 
Varieties (V) 19.33%*  8,.23** 9.97%  27.24"* 28.91%" 42.83%"  48.23** 13.62 
Planting Pattern (P) | 18.45** 6.00" —-7.65* = 16.63** 6.74 6.57% 9.52* 8.36% 
VxP 2.958 022". G51" - 26.53% O27" 0.16%S 0.0885 9,238 
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Figure 2. Changes in forage yield to 
different varieties and planting patterns. 
Note. Pl, P2 and P3 represent different 
planting patterns. 
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Figure 3. Changes in dry matter in response 
to different varieties and planting patterns. 
Note. Pl, P2 and P3 represent different 
planting patterns. 


Dry matter production of sorghum was 
significantly affected by cultivars and 
planting patterns (Figure. 3). Among 
cultivars, the cultivar Sorghum-2011 
produced highest dry matter of 19.25t ha’! 
followed by Jawar 2002, while cultivar JS- 
2002 (17.5t ha!) produced less dry matter. 
In case of planting pattern, the highest yield 
of dry matter (18.53 t ha‘') was recorded in 
P3 (40 x 30 cm). The minimum dry matter 
yield (18.0 t ha') was observed in Pi 
followed by Po. 
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Figure 4. Changes in dry percent in 
response to different varieties and planting 
patterns. Note. Pl, P2 and P3 represent 
different planting patterns. 
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Figure 5. Changes in crude protein (%) in 
response to different varieties and planting 
patterns. Note. Pl, P2 and P3 represent 
different planting patterns. 


The data regarding dry matter (%) of 
three sorghum cultivars affected by planting 
patterns is presented in figure 4. The 
maximum dry matter (%) was recorded in 
cultivar Jawar-2002 (34.13%), followed by 
Sorghum-2011, while minimum dry matter 
percentage (30.38 %) was observed in JS- 
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2002. Among different planting patterns, dry 
matter percentage was found significant, and 
it ranged from 32.05-36.19 %. The 
maximum dry matter percentage (36.19%) 
was recorded in the treatment in which 
sorghum was sown by using a P3 (40 x 30 
cm) planting pattern. The minimum dry 
matter percentage (32.05 %) was recorded in 
the treatment P; followed by P2. 


3.2 Qualitative Attributes 


The crude protein content of different 
cultivars as affected by various planting 
patterns is presented in figure 5. Maximum 
crude proteins (%) were observed in 
Sorghum-2011 (6.12%) followed by Jawar- 
2002 (6.03%), while minimum crude protein 
(%) was observed in JS-2002 (5.85%). In 
the case of planting patterns, crude protein 
percentages were found significant, and it 
ranged from 5.90-6.06% the highest crude 
protein (6.06%) was recorded in the 
treatment in P3 (40 x 30 cm) planting 
pattern. The minimum crude protein 
percentage (5.90%) was recorded in the 
treatment P;(60 x 20 cm) followed by P2. 


Analysis of variance indicates that the 
percentage of crude fiber was considerably 
diverse in all varieties of sorghum (Figure 
6). Maximum crude fiber (32.12%) was 
observed in Sorghum-2011, and it was 
followed by Jawar-2002 (30.38 %), while 
minimum crude fiber (28.33 %) was 
observed in JS-2002. The effect of planting 
pattern on crude fiber percentage was found 
significant, ranging from 28.9-32.8%. 


The maximum crude fiber percentage 
(32.88%) was recorded in treatment P3. The 
minimum crude fiber percentage (28.91%) 
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was recorded in the treatment P;, followed 
by Po. 
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Figure 6. Changes in crude fiber content 
(%) in response to different varieties and 
planting patterns. Note. Pl, P2 and P3 
represent different planting patterns. 
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Figure 7. Changes in crude protein (%) in 
response to different varieties and planting 
patterns. Note. Pl, P2 and P3 represent 
different planting patterns. 


The data regarding ash percentage shows 
significant variation among sorghum 
cultivars (Figure 7). Ash percentage was 
maximum (8.73%) in cultivar Sorghum- 
2011 which was followed by Jawar-2002 
(8.43%) and were minimum (7.98 %) in JS- 
2002. Among planting patterns, ash contents 
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were found significant and ranged from 
8.10-8.64 %. The maximum ash percentage 
(8.64 %) was recorded in P3(40 x 30 cm) 
planting pattern. The minimum ash 
percentage (8.10 %) was recorded in the 
treatment P; (60 x 20 cm) planting patterns 
followed by Po. 


4. DISCUSSION 


Selection of cultivar and planting pattern 
is one of the most important method to get 
higher yield and quality of sorghum crop. 
Therefore, we conducted a field experiment 
to evaluate the best responsive planting 
pattern and cultivar to improve sorghum 
yield and quality under _ sustainable 
agriculture. 

Our results showed that growth and yield 
parameters including planting density (m~”), 
plant height (cm), leaves plant', leaf area 
plant, stem diameter (cm), weight plant! 
(g) fresh weight plant! (g) and dry weight 
plant’ (g) were significantly increased in the 
planting pattern of P3 (40 x 30 cm) under 
sorghum 2011 cultivar. These increments 
might be due to differentiation in the genetic 
makeup of cultivars and the adaptability of 
these varieties to different environmental 
conditions (Blum, 2004). Yousif et al. 
(2012) also found that different sorghum 
cultivars vary in plant height. Similarly 
Zulfigar et al.(2009) documented that the 
number of leaves per plant were influenced 
by different planting pattern and cultivars. 
Present results contradict Miranda et 
al.(2013), who found non-significant 
dissimilarity in leaves number of various 
sorghum varieties. These contradictory 
results might have been due to differences in 
environmental conditions and the genetic 
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potential of the varieties. An increase in 
weight per plant of Sorgum-2011 was due to 
greater plant height, leaf area, and stem 
diameter. Similarly, Gondal et al.(2017) 
showed significant differences among 
sorghum cultivars regarding weight per 
plant. Furthermore, Ayub et al.(2010) found 
non-significant differences in plant dry 
weight among six sorghum varieties, while 
Afzal et al. (2012) compared different 
sorghum cultivars. They found a significant 
difference in weight per plant of different 
forage sorghum cultivars. 


Protein contents are a major parameter 
affecting forage crops’ nutritional value and 
quality. It’s the mixture of true protein and 
non-protein nitrogen and the fodder with 
high protein contents is considered a good 
quality fodder. While crude fiber percentage 
is one of the most important parameters that 
influences forage crops' quality. Fodders 
with low crude fiber contents are considered 
good quality forage because low fiber 
contents increase digestibility and 
palatability and improve intake. In addition, 
ash contents are described as the mineral 
contents in dry matter feedstuff and are 
mostly measured on a % basis. In the 
present study, protein contents, crude fiber 
percentage and ash content were improved 
in the planting pattern of P3 (40 x 30 cm) 
under sorghum 2011 cultivar. These 
significant differences may be due to the 
difference of growth stage at harvest. 
Another reason for these increments might 
be the more nutrient availability to plant in 
P3 patterns. Panwar et al. (2000) also 
showed significant differentiation in crude 
fiber contents among varieties of sorghum. 
Ayub et al. (2010) studied two sorghum 
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cultivars viz. JS -263 and Hegari for crude 
fiber percent and ash percentage revealed 
Hegari produced higher protein and ash 
contents than JS-263. 


5. CONCLUSION 


Based on our results, we observed that 
among different cultivars, sorghum 2011 
resulted in higher yield, yield components 
and quality traits as compared to other 
cultivars. Similarly, Sorghum-2011 also 
performed better at the Planting pattern of P3 
(40 x 30 cm) compared to other planting 
patterns. Therefore, we concluded that the 
cultivar Sorghum-2011 performed better 
than the other two cultivars because it 
produced a higher yield with quality at 
Planting pattern P3 (40 x 30 cm). It can be 
recommended for cultivation under 
Faisalabad conditions. 
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